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PREFACE 


This  report  is  a  product  of  the  Laboratory  Simulation  of  Spectral  and 
Directional  Spectral  Waves  Work  Unit,  Coastal  Flooding  and  Storm  Protection 
Program,  Civil  Works  Research  and  Development,  at  the  US  Army  Engineer  Water¬ 
ways  Experiment  Station's  (WES)  Coastal  Engineering  Research  Center  (CERC). 
Testing  was  conducted  from  February  to  April  1986,  and  data  reduction  and 
report  preparation  were  completed  in  April. 

This  report  was  prepared  by  Mr.  Michael  J.  Briggs,  Research  Hydraulic 
Engineer,  and  Ms.  Mary  L.  Hampton,  Civil  Engineering  Technician,  under  direct 
supervision  of  Mr.  Douglas  G.  Outlaw,  Chief,  Wave  Processes  Branch,  and  under 
general  supervision  of  Mr.  C.  E.  Chatham,  Chief,  Wave  Dynamics  Division, 

Mr.  Charles  C.  Calhoun,  Jr.,  Assistant  Chief,  and  Dr.  James  R.  Houston,  Chief, 
CERC.  Messrs.  John  H.  Lockhart,  Jr.,  and  John  G.  Housley,  Office,  Chief  of 
Engineers,  were  Technical  Monitors  for  the  Coastal  Flooding  and  Storm  Protec¬ 
tion  Program. 

Numerous  Individuals  contributed  to  successful  completion  of  this  proj¬ 
ect.  In  the  WES  Instrumentation  Services  Division  (ISD) ,  Mr.  Barry  W. 
McCleave,  Electronics  Engineer,  wrote  and  developed  much  of  the  wave  genera¬ 
tion  and  measurement  software  and  Electronic  Technicians,  Messrs.  David  A. 
Dailey  and  Lonnie  L.  Friar,  operated  and  maintained  the  directional  spectral 
wave  generator,  wave  gages,  and  associated  electronics.  Also  in  ISD, 

Messrs.  Homer  C.  Greer  III  and  Selwyn  W.  Guy  provided  trouble-shooting  ser¬ 
vices  when  electronic  problems  occurred.  In  CERC’s  Coastal  Processes  Branch, 
Mr.  Norman  W.  Scheffner,  Research  Hydraulic  Engineer,  developed  the  nonlinear 
cnoldal  waveform  generation  software  and  assisted  with  basin  bathymetry  mea¬ 
surements.  In  the  Wave  Processes  Branch,  Mr.  Kent  A.  Turner,  Computer  Pro¬ 
grammer  Analyst,  was  Invaluable  for  his  day-to-day  assistance  with  packages 
for  wave  analyses.  Mr.  Larry  A.  Barnes,  Civil  Engineering  Technician,  pro¬ 
vided  expertise  in  detailed  design  of  wave  absorbers  and  gage  support  frames, 
interfacing  with  the  WES  shops,  and  assisting  with  raport  preparation.  Con¬ 
tract  students,  Ms.  Christie  T.  Sanders  and  Mr.  Larry  D.  Davis,  generated  con¬ 
trol  signals  and  helped  with  data  reduction  and  report  preparation. 

Mrs.  Janie  G.  Daughtry  and  Mrs.  Dorothy  L.  Staer,  Branch  Secretaries,  prepared 
the  final  draft  report.  Ms.  Jamie  W.  Leach  of  the  WES  Information  Products 
Division  edited  this  report. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


_ Multiply _  _ _ 

feet  0.3048 

gallons  (US  liquid)  3.785412 

horsepower  (550  foot-pounds  745.6999 

(force)  per  second) 

Inches  2 . 54 

pounds  (mass)  0.4535924 

square  feet  0.09290304 


To  Obtain 
metres 

cubic  decimetres 
watts 

centimetres 
kilograms 
square  metres 


DIRECTIONAL  SPECTRAL  WAVE  GENERATOR  BASIN  RESPONSE 
TO  MONOCHROMATIC  WAVES 


PART  I :  INTRODUCTION 


1.  The  directional  spectral  wave  generator  (DSWG)  of  the  US  Army  Engi¬ 
neer  Waterways  Experiment  Station  (WES)  Coastal  Engineering  Research  Center 
(CERC)  is  a  unique  resource  for  study  of  natural  sea  states  In  a  laboratory 
environment.  It  is  the  only  one  of  its  kind  in  the  United  States  dedicated 
exclusively  to  simulation  and  analysis  of  naturally  occurring  short-crested 
waves  in  a  coastal  environment  and  wave  processes  associated  with  these  waves. 
Some  unique  features  are  size,  modular  design,  portability,  method  of  paddle 
connection  and  displacement,  and  electric  motor  power.  Basic  monochromatic 
waveforms  may  be  either  linear  sinusoidal  or  nonlinear  higher  order  cnoidal. 

In  addition  to  monochromatic  waves,  unidirectional  and  directional  sea  and 
swell  components  representative  ot  different  measured  and  empirical  wave  spec¬ 
tra  can  be  simulated.  Also,  wave  groups,  wave  transients,  and  other  wave 
forms  such  as  explosion  waves  can  be  simulated.  The  inclusion  of  directional 
spreading  in  design  calculations  for  offshore  and  coastal  structures  may  well 
result  in  less  costly  design.  Coastal  structures  respond  differently  to 
three-dimensional  wave  environments  than  to  two-dimensional  waves.  Three- 
dimensional  waves  create  movements  and  forces  transverse  to  the  main  wave  di¬ 
rection  that  affect  loads  (Danish  Hydraulic  Institute  (DHI)  1985).  Thus,  sig¬ 
nificantly  improved  physical  model  data  can  be  obtained  in  site-specific  and 
research  studies  using  the  DSWG  facility  and  naturally  occurring  sea  states. 

?.  In  late  1985,  the  DSWG  basin  was  relocated  from  its  original  posi¬ 
tion  (used  for  acceptance  testing  from  contractors)  for  more  efficient  utili¬ 
zation  of  model  space  in  the  movable-bed  test  facility  and  easier  access  to 
the  computer  control  room.  In  the  process,  a  new  wave  absorption  system  was 
designed,  built,  tested,  and  installed  along  the  basin  perimeter.  New  wave 
gage  support  frames  also  were  designed  and  built.  These  frames  allow  multiple 
gages  to  be  installed  in  a  linear  array  with  less  interference  from  support 
legs.  A  centrally  located  power  and  gage  connection  box  allows  greater  flexi¬ 
bility  for  relocating  the  DSWG  and  wave  gages.  Two  catwalks  located  parallel 
to  the  DSWG  permit  easy  overhead  photography  and  cinematography. 
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3.  Monochromatic  performance  tests  were  conducted  from  February  to 
April  1986.  A  series  of  111  monochromatic  sinusoidal  waves  was  generated  to 
determine  the  relationship  between  DSWG  control  signal  and  the  measured  waves 
at  nine  different  locations  in  the  wave  basin.  The  series  consisted  of  com¬ 
binations  of  five  wave  periods,  heights  (i.e.  generator  strokes),  and  direc¬ 
tions.  The  wave  periods  ranged  from  0.75  to  3  sec,  the  generator  strokes  from 
1  to  11  in.,*  and  the  directions  from  0  to  60  deg.  In  addition,  wave  heights 
just  prior  to  breaking  were  determined  for  the  five  wave  periods  by  varying 
control  signal  stroke  amplitudes.  Effects  of  wave  nonlinearity  also  were 
investigated  by  comparing  linear  sinusoidal  results  with  cnoidal  waveforms 
generated  using  second-order  generator  control  theory.  Thus,  the  purpose  of 
these  tests  was  to  verify  theoretical  predictions  and  correct  them  as  neces¬ 
sary,  to  ensure  repeatability,  to  note  any  differences  within  the  basin  due  to 
location  and  time,  and  to  study  the  effects  of  nonlinearity. 

4.  In  Part  II,  descriptions  of  the  DSWG  wave  basin  and  wave  generating, 
wave  measurement,  and  data  acquisition  and  control  systems  are  given. 

Part  III  describes  wavemaker  theory,  including  two-  and  three-dimensional 
height  transfer  functions,  maximum  prebreaking  wave  heights,  constancy  of  mean 
wave  height,  snake  principle,  spurious  wave  period  limits,  wave  nonlinearity 
parameters,  control  signal  generation,  test  design  and  averaging  criteria, 
harmonic  analysis,  and  wave  analysis  for  height,  period,  and  direction. 

Part  IV  documents  test  procedures  used  in  the  five  project  phases  including 
theoretical  predictions,  control  signal  generation,  wave  gage  calibrations, 
wave  generation  and  measurement,  and  analysis.  Test  results  are  presented  in 
Part  V,  and  effects  of  wave  maker  period,  stroke,  and  angle  of  propagation  on 
wave  profiles  are  shown.  Results  of  comparisons  between  predicted  and  mea¬ 
sured  values  for  wave  periods,  heights,  and  directions  are  given  in  tabular 
and  graphical  form.  Harmonic  analysis  of  measured  data  shows  the  total  per¬ 
cent  variance  contained  in  the  first  harmonic  components.  Also,  maximum 
strokes  required  to  generate  prebreaking  waves  as  a  function  of  wave  period 
are  reported.  The  constancy  of  wave  heights  within  the  measurement  area  and 
along  wave  crests  is  described.  Also,  effects  of  nonlinearity  on  waveform 
and  measured  period  and  height  are  discussed.  Finally,  Part  VI  contains  a 
summary  of  results  and  recommendations  for  future  research  and  improvements. 

*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  5. 
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PART  II:  GENERAL  DESCRIPTION  OF  FACILITIES 


5.  In  this  part,  descriptions  of  the  DSWG  wave  basin,  wave  generator, 
wave  measurement,  and  data  acquisition  and  control  systems  are  given. 

DSWG  Wave  Basin 


6.  The  DSWG  wave  basin  is  housed  in.CERC's  movable-bed  test  facility 
which  contains  approximately  70,000  sq  ft  of  laboratory  floor  space.  Figure  1 
is  a  schematic  illustrating  the  relative  position  of  the  basin  to  three  other 
basins  and  the  office/computer  complex. 


Figure  1.  Location  of  DSWG  basin  within  movable-bed  test  facility 
Physical  dimensions 

7.  The  DSWG  basin  is  96  ft  long  by  114  ft  wide  and  can  accommodate 
water  depths  to  2  ft  within  its  concrete  block  walls  (Photo  1).  The  basin  is 
within  a  larger  irregular  area  which  measures  138.3  ft  long  by  120.8  ft  wide 
(Figure  2).  Rather  than  erect  a  wall  parallel  to  the  DSWG  to  close  off  this 
area,  it  was  decided  to  incorporate  the  irregular  geometry  to  aid  in  wave 
absorption  and  energy  dissipation  behind  the  wave  absorbers.  The  basin  holds 
approximately  112,000  gal  of  water  when  filled  to  the  1-ft  level. 
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ABSORBER 


ROOF  SUPPORT  POSTS 


Figure  2.  DSWG  basin  and  coordinate  system 
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ABSORBER 


Basin  bathymetry 

8.  The  bathymetry  of  the  basin  was  measured  at  5-ft  intervals  using  a 
level  and  rod  in  the  70-  by  90-ft  area  in  front  of  the  DSWG.  Figures  3a 

and  3b  are  three-dimensional  surface  and  contour  plots  of  the  DSWG  basin, 
respectively,  at  interpolated  2. 5-ft  grid  points  on  a  90-  by  90-ft  grid.  A 
0.8  smoothing  factor  was  used  for  both.  Because  of  the  loss  of  resolution  due 
to  the  smoothing  factor,  both  surface  and  contour  plots  are  only  intended  to 
illustrate  relative  trends  and  not  absolute  values.  On  the  three-dimensional 
surface  plot,  the  origin  of  the  DSWG  coordinate  system  is  in  the  lower  right 
corner  (i.e.  DSWG  Paddle  1)  and  extends  along  the  bottom  to  the  center  of  the 
page  (i.e.  DSWG  Paddle  6U) .  On  the  contour  plot,  the  DSWG  extends  along  the 
bottom  of  the  plot.  The  maximum  variation  in  the  basin  bathymetry  is  1.08  in. 
The  highest  point  is  +0.51  in.  at  coordinates  X  =  35  ft,  Y  =  5  ft.  The  lowest 
spot,  measuring  -0.57  in.,  is  at  coordinates  X  =  70  ft,  Y  =  80  ft.  Standard 
deviation  for  the  data  set  of  285  measurements  is  0.18  in.  Appendix  A  con¬ 
tains  a  tabular  listing  of  the  measured  values. 

Wave  absorption  system 

9.  The  DSWG  basin  has  a  unique  beach/wave  absorption  system.  Porta¬ 
ble,  metal  frames  with  slopes  of  approximately  37  deg  are  installed  along  the 
basin  perimeter,  and  frames  behind  the  DSWG  have  two  beach  faces  with  individ¬ 
ually  adjustable  slopes.  Wave  absorption  and  energy  dissipation  are  provided 
by  two  layers  of  2-in.  horsehair  sandwiched  between  two  layers  of  expanded 
metal.  Photos  2  and  3  show  the  wave  absorber  behind  the  DSWG  and  the  beach 
wave  absorbers  opposite  the  DSWG,  respectively.  In  flume  tests  of  16  selected 
cases  covering  wave  periods  from  0. 75  to  3  sec  and  wavemaker  strokes  of  1  to 

9  in.,  an  average  reflection  coetficient  of  12.25  percent  was  measured  for  the 
horsehair  and  expanded  metal  absorbers.  Also,  the  irregular  basin  shape  be¬ 
hind  wave  absorbers  at  the  beach  end  tends  to  scatter  wave  energy  outside  the 
main  wave  basin.  The  DSWG  basin  is  located  inside  an  insulated  steel  frame 
metal  shelter  which  has  four  center  roof  support  posts  (8  in.  diam)  spaced  at 
40-ft  intervals  behind  the  beach.  Reflection  and  diffraction  of  energy  from 
these  support  posts  are  insignificant  because  of  their  location. 

Observation  platforms 

10.  Two  overhead  catwalks  parallel  to  the  DSWG,  one  over  the  far  end 
and  one  bisecting  the  basin,  are  excellent  vantage  points  for  photography  and 
video  documentation. 
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Wave  Generating  System 


11.  The  DSWG  was  designed  and  built  by  MTS  Systems  Corporation  of 
Minneapolis,  Minnesota.  It  is  the  only  one  of  its  kind  in  the  United  States 
and  measures  90  ft  long,  4  ft  2  in.  wide,  and  4  ft  2-7/16  in.  high  (Photo  4). 
It  consists  of  60  paddles  in  4  modules  (i.e.  15  paddles  per  module),  each 
2.5  ft  high  and  1.5  ft  wide.  The  paddles  operate  in  a  piston  (translational) 
mode  (more  efficient  and  representative  of  shallow-water  environments)  and  are 
driven  at  each  of  the  61  joints.  This  produces  a  cleaner,  more  continuous 
waveform  with  less  cross-wave  generation  for  wave  periods  greater  than  approx¬ 
imately  0.70  sec  than  does  driving  individual  paddles  (Sand  1979).  Wave  peri¬ 
ods  are  typically  above  0.6  sec.  The  range  of  strokes  is  ±6  in.  corresponding 
to  a  ±10  volt  input  signal.  Directional  waves  are  generated  using  the  "snake 
principle"  (Part  III) .  Offset  angles  between  paddles  can  be  continuously  var¬ 
ied  within  the  range  of  0  to  180  deg  to  produce  directional  waves  at  angles 
approaching  ±90  deg  for  most  wave  periods. 

12.  A  unique  feature  of  this  generator  is  portability.  The  modular 
design  (Figure  4)  allows  it  to  be  moved  to  different  locations  within  the 
movable-bed  test  facility  for  different  projects.  Photo  5  shows  the  front  and 
back  of  an  individual  DSWG  module.  Thus,  it  is  a  "wetback"  design  with  no 
bottom  or  end  seals.  Figure  5  is  a  schematic  of  the  two  drive  plates,  fixed 


DRIVE 
assembly/ 
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Figure  4.  Modular  design  of  DSWG 
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and  articulating  hinges,  and  flexible  plastic  plate  seals  which  slide  in 
guides  between  each  paddle  to  provide  continuity.  Four  lifting  points,  two  in 
front  and  two  in  back  of  each  7,250-lb  module,  are  provided  for  transporting 
the  DSWG.  Six  adjustable  mounting  pads  are  provided  for  leveling  each  module. 
Also,  a  catwalk  and  cable  tray  for  interconnecting  cabling  and  conduit  are 
provided. 

13.  The  DSWG  is  an  electronically  controlled,  electromechanical  system. 
Photo  6  is  a  close-up  of  the  individual  0.75-hp  electric  direct  current  (DC) 
motors  which  drive  each  paddle  joint.  Rotary  motion  of  these  motors  is  con¬ 
verted  to  a  forward/backward  motion  through  pulley  and  drive  assemblies.  Also 
contained  on  each  module  are  four  Hoffman  enclosures:  (a)  power  disconnect 
enclosure  for  the  power  breaker  switch  and  main  power  fuses  for  460-v,  60-Hz, 
3-phase  at  22-amp  power  (Photo  7) ;  (b)  power  distribution  enclosure  for  main 
power  with  four  transformers,  power  distribution  terminal  boards,  and  one 
convenience  outlet  (Photo  8) ;  (c)  pulse  width  modulated  (PWM)  servo  motor  con¬ 
trollers  enclosure  containing  three  PWM  controllers,  power  supply,  control¬ 
lers,  and  cooling  fans  for  power  and  signal  conditioning  (i.e.  Getty  ampli¬ 
fiers,  Photos  9  and  10);  and  (d)  Temposonics  transducer  modules  enclosure  for 
control  and  feedback  transducers  for  monitoring  paddle  position  (Photo  11). 
Also  located  on  each  module  are  three  400-ft-long  power  and  control  cables  and 
take-up  reels  for  electrical  Interface  with  the  system  control  console.  Cable 
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reel  1  controls  drive  disable  and  power  on/off.  Cable  reel  2  carries  the 
drive  command.  The  displacement  feedback  is  controlled  by  cable  reel  3. 

14.  A  wave  generator  control  console,  also  built  by  MTS,  is  located  In 
the  computer  control  room.  Photo  12  shows  front  and  rear  views  of  this  con¬ 
sole.  It  Includes  both  digital  and  analog  circuits  and  provides  digital  wave 
board  control  signals  from  the  VAX  11/750  computer  for  input  to  61  Preston  D/A 
converters  using  an  IEEE  488  interface  for  each  of  61  D/A  driven  servo  mod¬ 
ules.  This  console  weighs  875  lb  and  Is  5  ft  wide,  3.5  ft  deep,  and  6  ft 
high.  Bay  1  contains  a  calibration/test  Indicator,  three  MTS  450  servo/ 
controller  chassis  with  10  servo  controller  modules  (which  provide  drive  com¬ 
mand,  drive  enable/disable,  and  displacement  readout  signals)  each,  and  two 
connector  panels  for  controlling  modules  1  and  2.  Bay  2  is  similar  to  Bay  1 
with  electronics  for  controlling  modules  3  and  4.  An  additional  servo  con¬ 
troller  module  is  provided  for  a  total  of  31  in  Bay  2  (i.e.  61  total  for  sys¬ 
tem).  Also,  a  power  distribution  box  for  two  26-v  DC  power  supplies,  two  DC 
voltmeters,  one  input  connector  plug,  and  six  output  connector  plugs  are  in¬ 
cluded.  Finally,  Bay  3  contains  one  Preston  D/A  controller,  12  ft  of  IEEE  488 
cable,  four  D/A  converter  chassis  with  61  Preston  D/A  converter  modules,  the 
system  control  panel,  and  the  power  entry  panel  for  120-v,  60-Hz,  1-phase, 
22-amp  power.  Necessary  connections,  interface  cabling,  backplane  wiring, 
cooling  fans,  and  power  plug  mold  strips  are  provided  as  required.  A  function 
generator  is  provided  for  calibration  checking  of  the  system.  Photo  13  illus¬ 
trates  the  module  power  "on"  and  "off,"  emergency  stop,  master  reset,  and  pro¬ 
gram  "run"  and  "stop"  buttons  on  the  system  control  panel.  Also,  stroke  and 
displacement  error  limits  and  operating  range  adjustments  including  program 
span,  displacement  set  point,  error  gain  and  rate,  drive  zero,  and  error  zero 
are  provided  by  a  servo  module  panel  for  each  drive  (Photo  14).  Interested 
readers  should  refer  to  the  paper  by  Outlaw  (1984)  and  the  MTS  927.57  Instruc¬ 
tion  Manual  (MTS  System  Corporation  1984)  for  more  details  of  this  system. 

WaVe  Measurement  System 


Measurement  area 

15.  A  20-ft-sq  measurement  area,  located  parallel  to  and  10  ft  from  the 
DSWG  (symmetric  about  its  center  line),  was  selected  for  these  experiments. 
Nine  wave  gages  were  Installed  with  10-ft  spaclngs.  For  waves  generated  by 
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the  snake  principle,  the  best  measurement  area  is  as  near  the  generator  as 
possible,  but  a  minimum  distance  of  3  to  10  water  depths  in  front  of  it 
(Goring  and  Raichlen  1980)  to  avoid  diffraction  and  finite  length  effects. 
Also,  to  minimize  effects  of  reflected  energy,  the  measurement  area  should  be 
as  close  to  the  generator  as  possible  to  maximize  the  time  for  measurement  be¬ 
fore  reflected  waves  return.  As  wave  direction  increases,  usable  measurement 
area  decreases  because  limiting  wave  rays  leaving  the  ends  of  the  DSWG  are  at 
such  a  steep  angle  that  they  cut  through  the  measurement  area.  Figure  6  shows 
the  measurement  area  and  the  maximum  wave  direction  angle  of  49  deg  which 
could  be  generated  and  still  permit  measurement  by  all  gages.  For  the  60-deg 
wave  direction,  only  Gages  1,  2,  3,  and  6  could  be  used  for  measurements. 

Wave  gages 

16.  The  wave  gages  are  parallel  wire  resistance-type  sensors.  They 
measure  the  change  in  water  surface  elevation  with  time  and  record  it  as  a 
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voltage  which  is  later  converted  to  a  wave  elevation  tine  series  by  analysis 
software.  Because  resistance-type  gages  are  conductivity  probes*  they  are 
very  sensitive  to  electronic  and  temperature  drift  and  cleanliness  of  the 
water  (l.e.  dirt,  oil,  etc.).  Two  sizes  of  wave  gages  were  used  in  this 
study.  "Snail"  gages  have  a  maximum  calibration  range  of  3.25  in.  Larger 
"Jordan"  gages  have  a  calibration  range  of  10  in.  The  small  gages  provided 
better  resolution  where  wave  heights  did  not  exceed  their  range. 

Gage  support  frames 

17.  Three  aluminum  wave  gage  support  frames,  22  ft  in  length,  were 
designed  and  installed  in  the  measurement  area.  The  frames  reduce  possible 
interference  with  the  wave  train  by  support  legs.  Photo  15  gives  different 
views  of  these  frames. 


Data  Acquisition  and  Control  System 


18.  An  automated  data  acquisition  and  control  system  (ADACS)  is  used  to 
generate  and  transmit  control  signals,  monitor  DSWG  feedback,  and  analyze  and 
store  measured  data.  The  heart  of  the  ADACS  for  the  DSWG  is  a  DEC  VAX  11/750 
computer.  Table  1  summarizes  computer  specifications  and  system  components. 
The  DSWG  control  signal  is  updated  20  times  per  second,  and  data  for  each  mea¬ 
surement  channel  are  sampled  at  a  50-Hz  sampling  frequency. 
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PART  III:  WAVEMAKER  THEORY 


19.  The  wavemaker  theory  dlscuseed  In  this  part  uiumi  linear  wave 
theory  for  plane  piston  wavemakers  In  shallow  vater,  and  standard  assumptions 
that  the  fluid  is  incompressible  and  irrotatlonal  are  made.  Descriptions  of 
wavemaker  theory  including  two-  and  three-dimensional  height  transfer  func¬ 
tions,  maximum  prebreaking  wave  heights,  constancy  of  mean  wave  height,  snake 
principle,  spurious  wave  period  limits,  wave  nonlinearity  parameters,  control 
signal  generation,  test  design  and  averaging  criteria,  harmonic  analysis,  and 
wave  analysis  for  height,  period,  and  direction  are  given. 

Wave  Height  Transfer  Function 

20.  Blesel  (1954)  originally  darived  the  two-dimensional  transfer  func¬ 
tion  F2  describing  the  ratio  of  generated  wave  height  H  to  wavemaker 
stroke  S  .*  By  equating  the  volume  of  water  displaced  by  a  piston  wavemaker 
to  crest  volume  in  the  waveform,  Galvin  (1964)  showed  that  the  asymptotic  ap¬ 
proximation  for  this  transfer  function  Is  given  by 

F2  -  |  -  kh  kh  <  ir/ 10  (1) 


where 

k  ■  wave  number 
h  »  water  depth 

21.  Madsen  (1974)  and  Dean  and  Dalrymple  (1984)  describe  the  procedure 
used  to  obtain  the  complete  wavemaker  theory  for  plane  waves  produced  by  a 
piston  wavemaker.  The  boundary  value  problem  for  the  wavemaker  Involves  solu¬ 
tion  of  the  governing  velocity  potential  equation  (i.e.  Laplace  equation)  with 
appropriate  boundary  conditions.  These  conditions  are  the  dynamic  and  kine¬ 
matic  free  surface,  kinematic  bottom,  and  radiation  boundary  conditions .  By 
combining  these  into  a  kinematic  wavemaker  boundary  condition  describing  the 
horlsontal  displacement  and  velocity  of  the  wavemaker,  neglecting  higher  order 


*  For  convenience,  symbols  and  abbreviations  are  listed  In  the  Notation 
(Appendix  G).  An  effort  has  been  made  to  conform  with  the  terminology 
recoaawnded  by  the  IAHR  (International  Association  for  Hydraulic  Research) 
Working  Group  on  Wave  Generation  and  Analysis  (1985). 
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terms  (i.e.  linearizing)  and  Invoking  orthogonality  conditions,  the  complete 
two-dimensional  transfer  function  F2  can  be  shown  to  be 


_  ,  .  H  2  cosh  (2kh  -  1) 
Vr;  ’  S  "  slnh  C2kh)  +  2kh 


(2) 


An  equivalent  form,  based  on  trigonometric  substitutions,  is  given  by  Sand 
(1979)  and  Sand  and  Lundgren  (1981)  as 


F2(f) 


H  _  2  slnh2  (kh) 

S  slnh  (kh)  cosh  <!kh)  +  kh 


(3) 


Equations  2  and  3  are  valid  over  all  ranges  of  nondlmenslonal  depth  kh  . 

22.  Standing  waves  which  the  wavemaker  may  create  due  to  mismatch  with 
the  water  particle  horizontal  velocity  are  assumed  to  decay  rapidly  with  dis¬ 
tance  from  the  wavemaker .  Thus,  measurements  are  usually  made  several  water 
depths  or  wavelengths  In  front  of  the  wavemaker. 


Directional  Effect  on  Wave  Height  Transfer  Function 

21.  The  three-disienslonal  wave  height  transfer  function  F^  Includes 
the  effects  of  wave  directionality  on  wave  height.  Using  Huygens  principle 
end  energy  flux  arguments.  Sand  (1979)  showed  how  F^  Is  an  extension  of  F^ 
for  multiple,  piston  type.  Infinitely  small  paddles.  To  first  order,  F^  Is 
given  by 


Fj(f ,0) 


r2(f) 

cos  6  S 


(4) 


where 

6  ■  wave  direction 
H(0)  -  corrected  wave  height 

Maximum  Ptcbreaklna  Wave  Heights 


24.  For  piston  wavemokera,  the  maximum  monochromatic  wave  height  which 
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can  be  generated  at  the  vavemaker  has  been  empirically  derived  by 


Ahrens*  as: 


H.  »  0.1L  tanh  (kh)  (5) 

D 

where  L  ■  linear,  shallow-water  wave  length.  Breaking  of  the  waves  at  the 
wavemaker  tends  to  occur  if  a  constant  larger  than  0.1  is  used  in  Equation  5. 
This  relationship  does  not  rule  out  the  possibility  of  waves  larger  than  this 
occurring  at  points  within  the  basin  as  a  result  of  shoaling,  etc.  It  only 
provides  that  they  do  not  break  at  the  wavemaker  because  of  the  discontinuity 
of  the  paddle  motion.  Equation  5  is  in  agreement  with  Naeser  (1981)  for  a 
practical  breaking  limit  H/L  >0.1  for  wave  periods  less  than  1.6  sec. 

Constancy  of  Mean  Wave  Height 

25.  The  mean  wave  height  variation  AH  is  a  measure  of  the  constancy 
of  wave  height  along  a  wave  front  within  the  measurement  area.  It  is  defined 
in  percent  as 


AH  -  100 


y  ("n  -  "> 

n-L  95 


whore 


H  -  wave  height  measured  at  wave  gage  n 
n 


H  -  average  wave  height  of  the  nine  gages  within  the  measurement  area 
Sand  (1979)  found  a  range  in  AH  In  his  experiments  of  AH  <  15  percent  for 
directions  0  s  10  deg  and  AH  <  25  percent  for  10  <  6  s  25  deg. 

26.  Possible  sources  of  wave  height  variation  in  wave  basins  are  the 
interaction  of  free  and  bound  higher  harmonics,  reflections,  end  diffraction, 
Benjamin-Fair  instability,  basin  assynsMtry,  variable  wavemaker  strokes,  vari¬ 
able  bottom  gap,  and  calibration  errors.  In  nature,  waves  are  characterized 
by  a  higher,  sharper  crest  and  flatter  trough  than  predicted  by  linear. 


*  Personal  Communication,  January  1986,  John  P.  Ahrens,  Research  Oceanogra¬ 
pher,  US  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  Miss. 
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first-order  wave  theory.  This  is  due  to  a  second-order  effect  of  a  bound  or 
locked  higher  harmonic  (second  harmonic  will  have  greatest  amplitude)  wave 
(BHW)  which  consists  of  sums  of  fundamental  frequencies.  If  first-order  wave 
theory  is  used  for  calculating  the  wavemaker  control  signal,  boundary  condi¬ 
tions  at  the  wavemaker  cannot  be  satisfied  for  the  second-order  BHW,  resulting 
in  the  creation  of  spurious  or  free  higher  harmonic  waves  (FHW) .  The  result 
in  a  wave  basin  is  a  continually  changing  wave  profile  as  a  function  of  posi¬ 
tion  due  to  alternating  cancellation  and  reinforcement  with  the  BHW  component. 
As  the  water-depth-to-wave-length  ratio  decreases,  a  sinusoidal  control  signal 
transforms  in  a  basin  into  a  nonlinear  waveform  by  binding  a  second  harmonic 
BHW  to  the  fundamental  wave.  In  the  process,  an  FHW  is  liberated  which  dis¬ 
torts  the  wave  profile  as  described  above  (Sand  and  Mansard  1985;  Briggs  1986). 

27.  Reflections  are  another  major  contributor  to  wave  height  variabil¬ 
ity  in  wave  basins.  A  reflection  coefficient  of  10  percent  results  in  a  wave 
height  variation  of  ilO  percent.  Reflections  are  greatest  for  waves  with  pe¬ 
riods  larger  than  2  sec  because  low  frequency  energy  is  less  easily  absorbed. 
Also,  during  reflection,  wave  energy  is  transferred  between  BHW  and  FHW  waves 
as  the  fundamental  wave  is  partly  absorbed,  partly  reflected,  and  partly  con¬ 
verted  to  FHW  waves  (Mansard,  Sand,  and  Funke  1985). 

28.  Energy  can  be  diffracted  from  the  ends  of  a  wavemaker  that  is  not 
sealed  at  the  ends.  These  diffracted  waves  affect  the  measurement  area  di¬ 
rectly  and  also  reflect  off  side  boundaries  into  the  Measurement  area. 

29.  According  to  Funke,*  low  frequency  standing  cross  waves  due  to  mean 
water  level  changes  in  a  shallow-water  basin  might  act  to  modulate  a  basic 
sine  wave  component  causing  shedding  of  energy  to  different  frequency  bands, 
producing  a  type  of  "Benjemin-Felr"  effect. 

30.  A  variable  bottom  gap  between  the  bottom  of  individual  paddles  in 
the  wavemaker  and  the  wave  basin  floor  will  produce  smaller  and  more  spatially 
variable  wave  height  than  theoretically  predicted.  Also,  if  transducers  con¬ 
trolling  the  amount  of  stroke  going  to  each  paddle  are  out  of  calibration, 
slight  variability  in  generated  wave  heights  could  result.  According  to 
Madsen  (1971,  1974)  leakage  and  turbulence  adjacent  to  the  wavemaker  produce 
variability  in  wave  height. 

*  Personal  Communication,  29  May  1986,  E.  R.  Funke,  Senior  Research  Officer, 
Hydraulics  Laboratory,  Rational  Research  Council  of  Canada,  Ottawa,  Ontario, 
Canada. 
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31.  Finally,  assymmetry  of  Che  wave  basin  can  produce  complicated  re¬ 
flection  patterns  as  the  beach  response  Is  not  uniform.  Also,  calibration 
tolerances  of  the  wave  gages  will  affect  measured  values.  Thus,  it  is  normal 
to  expect  a  certain  amount  of  variation  in  measured  wave  heights  due  to  a  com¬ 
bination  of  one  or  more  of  these  causes. 

Snake  Principle 

32.  The  DSWG  utilizes  the  snake  principle  to  generate  waves  with  direc¬ 
tions  approaching  190  deg.  This  terminology  evolved  because  the  wavemaker 
looks  more  or  less  like  an  undulating  snake  as  the  waveform  progresses  along 
its  length.  For  monochromatic  waves,  energy  is  generated  at  a  discrete  fre¬ 
quency  and  direction  using  a  constant  phase  lag  or  offset  angle  between  pad¬ 
dles.  Several  authors  have  described  the  procedure  Including  Naeser  (undated) 
and  Offshore  Technology  Corporation  (1984).  The  angle  of  wave  propagation  or 
wave  direction  0  is  determined  by 


sin  0  *  ^  (7) 

where  Y  *  distance  along  the  wavemaker  corresponding  to  one  period  or  cycle 
of  360  deg  (Figure  7).  The  distance  Y  is  related  to  the  paddle  width  B  by 

Y  -  N^B  (8) 

where  ■  number  of  paddles  required  to  make  one  cycle  and  is  a  function  of 
the  offset  phase  lag  6 


u  .  360 

Nt  T 


(9) 


Figure  8  illustrates  how  the  angle  of  wave  propagation  changes  as  the  sign  of 
the  phase  angle  changes. 

33.  Naeser  et  al.  (1981)  and  Sand  (1979)  note  that  the  generation  of 
"snake  waves"  is  restricted  for  periods  below  1  sec  because  of  the  need  to 
have  a  minimum  of  two  paddles  over  the  cycle  distance  Y  (i.e.  phase  lag  of 
180  deg).  Thus,  for  the  0. 75-sec  period,  it  is  not  possible  to  generate  waves 
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Figure  7.  Snake  principle  for  generating  waves 


WAVEMAKER 


Figure  8.  Direction  of  wave  propagation  relative  to 

phase  angle 


up  to  ±90  deg.  A  more  realistic  maximum  wave  direction  angle  is  approximately 
70  deg  for  an  offset  phase  of  179.1  deg.  For  an  integer  number  of  paddles  as 
used  in  these  tests,  the  maximum  wave  angle  declines  to  38.8  deg  for  this  wave 
period.  Table  2  lists  the  values  for  the  offset  and  direction  angles  for  ex¬ 
act  and  integer  number  of  paddles  over  a  wavelength.  The  maximum  wave  direc¬ 
tion  angles  for  an  integer  number  of  paddles  were  used  in  this  study.  The 
DSWG  can  generate  waves  within  a  continuous  range  of  angles  as  long  as  two  or 
more  paddles  are  used. 


Spurious  Wave  Period  Limits 

34.  As  a  consequence  of  the  finite  width  of  the  wavemaker  paddles, 
small  undesired  disturbances  known  as  spurious  waves,  distinct  from  the  spuri¬ 
ous  waves  discussed  in  paragraph  26,  are  generated  in  addition  to  the  main 
waves.  The  number  of  these  waves  is  fixed  by  paddle  width  and  wavelength,  and 
their  frequency  is  the  same  as  that  of  the  main  wave.  Amplitude  and  direc¬ 
tions,  however,  change  for  each  spurious  wave  component  p  .  Spurious  wave 
amplitudes,  small  relative  to  the  main  wave,  are  attenuated  even  more  as  the 
number  of  spurious  wave  components  increases.  The  DSWG  is  driven  at  the 
joints  between  the  paddles  rather  than  at  the  center  of  individual  paddles 
(Figure  5) .  Although  this  does  not  change  the  number  of  spurious  waves  which 
may  be  generated.  Sand  (1979)  and  Sand  and  Lundgren  (1981)  showed  that  it  does 
greatly  attenuate  the  amplitude  of  waves  which  might  be  produced. 

35.  According  to  Naeser  et  al.  (1981),  the  minimum  period  T  neces- 

min 

sary  to  prevent  spurious  waves  is  defined  by 


T  ,  -  0.53  -J  /Z  +  sin  0  (10) 

min  * 

where  0  =  angle  between  the  wave  crest  and  the  wavemaker.  Figure  9  illus¬ 
trates  this  relationship.  Table  3  is  based  on  Sand’s  results  and  lists  the 
number  of  spurious  waves  generated  for  wave  periods  less  than  1.0  sec  with  the 
range  of  wave  directions  given.  Because  waves  with  periods  less  than  0.75  sec 
were  not  considered  in  this  study,  spurious  wave  generation  did  not  pose  a 
serious  problem. 
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Figure  9.  Minimum  spurious  wave  periods 

36.  The  wave  height  of  the  pth  spurious  wave  component  H  is 

P 


H  -  H(0) 
P 


8ln  fc) 

’(’  - 


p  *  1,2... 


where  H(0)  -  wave  height  of  the  main  wave  corrected  for  directional  effects 
(see  paragragh  23  on  directional  effects  on  wave  height  transfer  function). 

37.  Sand  also  showed  that  the  direction  0^  of  the  ptfl  spurious  wave 
is  given  by 


P  "  1,2,... 


Wave  Nonlinearity  Parameters 


38.  Certain  wave  height  and  period  combinations  produce  waves  which 
tend  to  deviate  from  the  linear  Airy  wave  profile  (i.e.  sinusoidal).  The 
waves  tend  to  become  more  nonlinear:  higher,  sharper  crests  and  broader, 
flatter  troughs.  Higher  harmonic  components  become  more  pronounced  in  the 
wave  profile  on  the  back  side  of  the  crest  and  in  the  troughs.  This  phenome¬ 
non  is  due  to  mismatch  in  wavemaker  motion  with  the  required  water  particle 
velocities.  As  the  waves  become  more  nonlinear,  it  is  necessary  for  the  wave- 
maker  to  move  forward  faster,  and  slower  coming  back.  Descriptions  of  the 
control  signal  generation  for  both  linear  and  nonlinear  cnoidal  control  theory 
waveforms  are  given  in  the  next  section. 

39.  Measures  of  wave  nonlinearity  are  the  wave  steepness  H/L  and 
Goda's  (1985)  nonlinearity  parameter  II  defined  as 

n  -  jr  coth3  (kh)  (13) 

Whereas  the  Ursell  number  is  only  appropriate  for  shallow-water  conditions, 
Goda's  parameter  can  be  used  for  all  water  depths  and  reduces  to  Ursell' s  for 
shallow  water.  It  appears  that  waves  with  values  of  II  less  than  0.1  are 
relatively  linear  with  the  nonlinearity  increasing  above  this  value.  Compari¬ 
sons  of  wave  period  and  height  with  this  parameter  are  made  in  Part  V. 


Control  Signal  Generation 


40.  The  primary  type  of  control  signal  was  a  linear  sinusoidal  wave. 

To  ascertain  effects  of  nonlinearity  on  higher  waver  periods,  nonlinear  cnoi¬ 
dal  control  signals  also  were  generated.  Orientation  of  the  coordinate  system 
is  the  same  as  in  Figure  2  with  the  origin  at  the  outer  edge  of  paddle  1,  the 
x-axls  extending  perpendicular  to  the  wavemaker,  and  the  y-axis  parallel  to 
the  wavemaker. 

Linear  sinusoidal  waveforms 

41.  For  the  sinusoidal  waveform,  the  governing  equation  for  control 

signal  to  each  of  the  61  pistons  at  location  y  and  time  t  is 


04) 


Sc(y,t) 


cos  (2irft  +  $  ) 


where  stroke  S  ,  frequency  f  ,  and  phase  <f>y  are  input  quantities.  The 
wave  height  at  a  particular  location  in  the  basin  is  a  function  of  the  three- 
dimensional  stroke-to-height  transfer  function  previously  discussed. 

42.  The  phase  consists  of  a  constant  portion  and  an  offset  portion 

which  determines  the  wave  angle  generated.  The  constant  portion  was  taken  as 
zero.  The  wave  is  propagated  along  the  wavemaker  at  the  proper  angle  by  off¬ 
setting  successive  paddles  by  the  offset  phase.  This  offset  portion  is  main¬ 
tained  for  the  entire  time  of  wave  generation.  The  phase  automatically  re¬ 
cycles  every  360  deg. 

43.  The  stroke  time-history  is  converted  to  a  voltage  time-history 
using  the  wavemaker  conversion  factor.  Since  20  v  *  12  in.  for  the  maximum 
possible  stroke  and  the  A/D  converter  has  a  resolution  of  4,096  digital  units, 
the  corresponding  factor  is  341.33  digital  units/in.  or  a  resolution  of 
0.00293  in. 

Nonlinear  cnoidal  waveform 

44.  In  the  generation  of  shallow-water  long  waves  with  constant  shape, 
Goring  and  Raichlen  (1980)  noted  that  it  is  important  to  consider  nonlinear 
aspects.  His  method  eliminates  the  trailing  waves  ("oscillatory  tail")  usu¬ 
ally  present.  The  velocity  of  the  paddle  as  it  moves  is  matched  with  the 
depth-averaged  water  particle  velocity.  The  governing  equation  for  the  paddle 
displacement  time-history  S(t)  is 


5(t) 


-  f 

kh  J 


f(w)  dw 


(15) 


where 

0  -  k(Ct  -  K) 

C  *  celerity 
f(w)  ■  function  of  ® 

w  -  dummy  variable  of  integration 

45.  A  cnoidal  waveform  time  series  with  360  Increments  is  generated 
based  on  Input  water  depth  and  desired  wave  height  and  wavelength.  Because 
the  cnoidal  waveform  is  different  from  the  sinusoidal,  period  and  wavelengths 
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do  not  match  exactly  with  those  predicted  by  linear  wave  theory.  For  unifor¬ 
mity,  the  same  wave  periods  were  used  for  the  cnoidal  wave  control  signals, 
and  the  wavelengths  were  allowed  to  change  as  necessary. 

46.  Each  of  the  360  elements  of  the  cnoidal  waveform  are  matched  to  the 
appropriate  phase  for  each  paddle  at  each  time  step.  The  total  phase 

<J)^  «  2irft  +  <t>  is  the  same  one  described  in  the  paragraphs  above.  The  first 
part  controls  the  frequency  increment  over  a  cycle,  and  the  second  part  con¬ 
sists  of  a  constant  and  offset  component.  Again,  the  constant  portion  is  set 

to  zero,  and  the  offset  is  a  function  of  the  wave  angle  desired. 

47.  Finally,  the  control  signals  for  all  61  pistons  are  converted  to 
the  corresponding  voltage  time  series  using  the  procedure  previously  described. 

Test  Design  and  Averaging  Criteria 

48.  In  defining  the  basin  response  to  different  wave  conditions,  it  is 

preferable  to  obtain  as  many  averages  of  the  measured  parameters  as  possible. 

All  waves  averaged  need  to  be  uniform  in  quality.  At  the  beginning  of  wave 

generation,  the  DSWG  imposes  a  10-sec  ramp  time  to  protect  system  components. 
Also,  a  minimum  number  of  waves  must  be  generated  to  ensure  that  enough  time 
has  elapsed  for  the  waves  to  travel  to  the  wave  gages.  As  the  wave  generation 
continues,  the  Benjamin-Feir  (Benjamin  and  Feir  1966)  instability  may  cause 
the  waves  to  oscillate  in  amplitude  with  time.  Also,  reflected  energy  could 
contaminate  the  measured  signal.  Although  flume  tests  over  the  range  of  peri¬ 
ods  used  indicated  reflection  coefficients  averaging  12.3  percent  for  the  new 
wave  absorber  beach,  any  reflected  energy  would  tend  to  degrade  the  incident 
waves.  Therefore,  test  cases  were  run  to  determine  the  maximum  number  of 
waves  which  could  be  averaged  for  each  period.  The  predicted  travel  times  for 
the  first  wave  to  go  to  the  back  gage  row  and  to  reflect  off  the  back  wall  and 
return  were  compared  with  measured  values. 

49.  Figures  lOa-lOe  illustrate  the  travel  paths  of  wave  rays  from  the 
left  and  right  ends  of  the  wavemaker  through  one  complete  reflection  for  each 
angle  of  wave  propagation,  0  to  60,  respectively.  The  travel  time  Tj  to  the 
back  gage  row  located  30  ft  from  the  wavemaker  is  given  by 


T 


1 


30 


C  cos  0 

g 


(16) 
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where  -  group  velocity, 
back  wall  (i.e.  134  ft  -  82 
back  wall  to  back  gage  row) 


The  analogous  travel  time  to  reflect  off  the 

ft  from  wavemaker  to  back  wall  plus  52  ft  from 
and  return  is 


.  134 

T2  "  C  cos  0 

g 


(17) 


Even  though  the  wave  ray  may  undergo  several  reflections  prior  to  returning  to 
the  back  gage  row  (i.e.  45  deg  and  60  deg),  the  total  distance  traveled  paral¬ 
lel  to  the  x-axls  is  the  same. 


Zero-Crossing  and  Harmonic  Analysis 


Zero-crossing 

50.  Two  types  of  wave  analysis  are  performed  on  the  equally  spaced  mea¬ 
sured  data:  zero-crossing  and  harmonic  analysis.  Zero-crossing  analysis 
calculates  average  and  significant  wave  heights  and  periods  and  their  standard 
deviations  using  a  windowing  technique  based  on  the  specified  period.  Each 
windowed  segment  of  the  data  is  searched  for  the  minimum  and  maximum  value, 
assuming  that  each  crest  must  follow  a  trough  and  the  calculated  period  is 
within  a  specified  tolerance  of  the  generated  period.  Figure  11  is  a  flow¬ 
chart  of  the  search  procedure  from  a  report  by  Turner  and  Durham  (1984).  Dur¬ 
ing  the  search  process,  the  raw,  unsealed.  Integer  voltages  are  used  because 
Integer  arithmetic  is  faster  and  requires  less  memory  than  floating  point  num¬ 
bers.  The  data  are  later  scaled  by  the  calculated  calibration  coefficients 
using  a  linear  or  quadratic  fitting  procedure.  This  type  of  analysis  is  very 
sensitive  to  choice  of  window  length  and  tolerance  on  the  wave  period 
specified. 

Harmonic 

51.  The  least-squares  harmonic  analysis  is  based  on  a  Legendre  method 
to  fit  a  series  of  four  components  (for  this  analysis)  to  the  measured  wave 
elevation  time  series  (Turner  and  Durham  1984).  The  purpose  of  this  analysis 
is  to  ascertain  the  "purity"  of  the  measured  waveform  by  calculating  the 
amount  of  total  energy  contained  in  the  fundamental  frequency  and  three  higher 
harmonic  components.  Ideally,  100  percent  of  the  variance  is  contained  in  the 
fundamental  frequency  component.  If  a  linear  control  signal  is  used  to 
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Figure  10.  (Sheet  2  of  3) 
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e.  Wave  direction  *  60  deg 
Figure  10.  (Sheet  3  of  3) 
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T 


generate  e  wave  which  Is  highly  nonlinear,  energy  will  be  shifted  to  harmonics 
of  the  fnndnaental  In  «n  attempt  to  produce  the  nonlinear  wave  profile. 

52.  The  Measured  water  surface  elevation  tine  aeries  na(t)  can  be 
assumed  to  be  periodic  In  fora  and  composed  of  a  true  component  n((t)  and  a 
noise  component  e(t)  given  by 

J 

n^Ct)  -  (t )  ♦  c(t)  “  *0  +  ^  *j  coe  (<**jt)  +  bj  aln  (w^t)  +  t(t)  (18) 

J-l 


where 


and  b. 


J  “  total  nuaber  of  coaponenta 
Fourier  coefficients 


angular  frequency  of  the  jth  coaponent 


The  noise  level  la  not  known  and  the  unknown  Fourier  coefficients  are  deter- 
alned  by  a  least-squares  procedure  which  minimizes  the  variance  E  such  that 


N  N 

E  -  £  c2(nAt)  -  £  kMt)  -  na(nAt)]2  (19) 

n**l  n«l 


where 

N  -  total  nuaber  of  data  saaples 

At  ■  tlae  Interval  (l.e.  0.02  sec)  between  consecutive  saaples 
To  minimize  the  variance,  calculate 


J  “  1 ,2, . . . , J 


(20) 


and 


3E 

*7 


0 


Equation  19  results  In  a  set  of  2J  simultaneous  equations  that  may  be  solved 
for  the  Sj  and  b^  coefficients.  It  may  be  rewritten  as 
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(21) 


lt(t)  -  Aq  +  Aj  cos  (oijt  +  ♦j) 

J-l 


where  the  amplitude  A^  and  the  phase  4^  of  the  j  component  are  deter¬ 
mined  froa  the  Fourier  coefficients 


1*1  +  bJ 


and  Aq  “  swan  water  surface  elevation.  Thus,  a  true  or  estimated  surface  el¬ 
evation  tiae  series  is  calculated  from  which  percentage  of  total  variance  for 
each  frequency  can  be  determined. 

Wave  Direction  Analysis 

53.  The  direction  of  wave  travel  is  calculated  by  computing  travel  time 
for  a  wave  front  to  travel  between  any  two  gages  in  the  array.  This  travel 
tiae  can  be  calculated  using  cross-correlation  techniques  or  measured  manu¬ 
ally.  In  this  study,  the  latter  approach  was  used. 

54.  Figure  12  illustrates  the  procedure  for  calculating  wave  direction. 
The  travel  tiae  Ac  for  a  wave  to  go  froa  Gage  2  to  Gage  3  is  measured  from 
the  respective  tiae  series  plots  based  on  phase  lag.  The  Incremental  wave¬ 
length  Al  then  is  given  by 

AL  -  CAt  (23) 


The  measured  wave  direction  3  then  is 

a 


a  .  -1  / Al\ 

0  -  sin  l  r-  ) 

a  \  fly  f 


where  Ay  -  distance  between  gages  parallel  to  the  waveaaker. 


34 


DIRECTION  OF  PROPAGAT 


Figure  12.  Calculation  of  wave  direction  angle 
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PART  IV:  TEST  PROCEDURES 


55.  Descriptions  of  the  test  procedures  are  contained  In  five  phases: 

a.  Theoretical  predictions. 

b.  Control  signal  generation. 

c.  Wave  gage  calibration. 

d.  Wave  generation  and  measurement. 

e.  Wave  analysis. 

Each  of  these  Is  discussed  in  the  paragraphs  which  follow. 


Theoretical  Predictions 


56.  In  this  Initial  phase,  certain  properties  based  on  linear,  depth- 
corrected,  Airy  wave  theory  were  calculated  for  each  of  the  111  monochromatic 
waves  desired.  These  properties  were  used  in  the  design  and  planning  of  tests 
and  generation  of  control  signals.  They  are  listed  in  Table  4  in  groups  ac¬ 
cording  to  their  functions. 

57.  The  computer  program  MONOSUMMARY  was  written  to  perform  these  cal¬ 
culations.  It  consists  of  a  main  driver  program  and  17  subroutines  described 
In  Table  5.  The  program  hierarchy  is  straightforward:  all  subroutines  are 
called  from  the  main  program  in  the  order  listed.  The  basic  input  parameters 
are  water  depth,  wave  period,  generator  stroke,  and  wave  direction.  Table  6 
lists  these  and  the  other  Input  variables  required.  An  example  of  command 
procedure  to  run  the  program  Is  contained  In  Appendix  B.  Descriptions  of  the 
output  variables  are  given  In  Table  /.  finally,  results  from  calculations  for 
each  of  the  111  test  cases  are  also  contained  In  Appendix  B. 


Control  signal  Generation 

58.  Either  sinusoidal  or  cnoldal  waveforms  can  be  created  In  this 
phase.  The  program  COMPONM  can  create  one  or  more  (maximum  of  256)  sinusoidal 
components  which  can  be  superimposed  to  obtain  a  desired  monochromatic  or 
spectral  sea  state.  Program  Inputs  Include  wave  amplitude,  period,  direction, 
and  phase  (see  Part  III).  Likewise,  program  C0Mf0NC4  can  create  multiple 
cnoldal  wave  components.  The  user  is  prompted  for  desired  wave  height,  wave¬ 
length  (period),  and  water  depth.  The  control  signal  can  be  started  at  the 
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origin  or  center  of  the  generator.  Each  paddle  is  advanced  in  phase  by  a  fre¬ 
quency  increment  corresponding  to  a  20-Hz  D/A  digitization  rate.  Thus,  after 
completion  of  a  cycle,  the  control  signal  automatically  recycles.  It  can 
continue  indefinitely  until  terminated  by  a  control  "C."  The  user  then  is 
prompted  to  close  the  file  and  tape  with  two  end-of-file  (EOF)  marks.  Appen¬ 
dix  Tables  Cl  and  C2  contain  listings  of  the  control  files  for  programs 
COMPONM  and  C0MP0NC4,  respectively.  Also,  Tables  Cl  and  C2  summarize  the 
operating  procedure  for  each  program,  respectively. 

Wave  Gage  Calibration 

59.  The  process  IDCAL  is  used  to  calibrate  wave  gages  and  ensure  that 
proper  gage  potentiometer  coefficients  are  used.  It  also  contains  descriptive 
information  for  documenting  test  output  and  data  files  for  archives.  Individ¬ 
ual  functions  are  organized  in  one  of  15  headers  which  permit  input  of  various 
parameters  including  number,  name,  coordinates,  and  relay  numbers  of  the 
gages,  potentiometer  and  rod  coefficients,  mean  depth  record,  etc.  Turner  and 
Durham  (1984),  Briggs,  Scheffner,  and  Hammock  (1985),  and  Hampton  (1986) 
describe  these  options  in  more  detail.  Prior  to  taking  data,  IDCAL  is  exe¬ 
cuted  to  place  correct  calibration  coefficients  in  a  generic  file  named 
"DUA1 : [DATA_STOjFOR018.DAT".  Whatever  coefficients  are  in  the  input  file 
specified  are  used  unless  the  gages  are  recalibrated  in  this  current  run. 

These  coefficients  then  are  used  in  the  analysis  to  convert  measured  voltages 
to  inches.  Appendix  Table  D1  describes  the  procedure  for  using  IDCAL  when 
calibration  coefficients  from  a  previous  run  are  used.  Table  8  lists  the  val¬ 
ues  used  in  header  2  for  each  of  the  five  wave  periods  tested. 

60.  At  the  beginning  and  end  of  each  day  of  testing,  wave  gages  were 
calibrated  according  to  the  procedure  described  in  Appendix  Table  D2.  It  is 
assumed  that  parameters  specified  in  the  other  headers  in  process  IDCAL  have 
been  previously  specified.  Wave  gages  are  calibrated  by  physically  moving 
gage  sensor  rods  through  a  series  of  11  steps  to  obtain  calibration  coeffi¬ 
cients  using  a  least-squares  linear  or  quadratic  fit.  This  averaging  tech¬ 
nique,  using  21  voltage  samples  per  gage,  minimizes  the  effects  of  slack  in 
gear  drives  and  hysteresis  in  the  sensors.  Appendix  Table  D3  lists  the  qua¬ 
dratic  fit,  maximum  deviation  calibration  coefficients  (in  units  of  feet  times 
10  ^)  for  each  of  the  small  gages  for  each  day  of  testing.  Also  included  are 
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the  average,  minimum,  and  maximum  values.  Appendix  Table  D4  contains  corre¬ 
sponding  values  for  Jordan  wave  gages  on  days  when  they  were  used. 

61.  The  constancy  of  the  1-ft  water  depth  was  verified  every  day  by 
taking  a  reading.  A  water  level  float  and  control  valve  were  installed  to 
maintain  the  depth  within  a  tolerance  of  ±0.001  ft  (0.012  in.)  or  0.10  percent 
after  8  March.  Prior  to  this  date,  the  water  level  was  manually  checked  and 
adjusted  each  day  to  within  a  tolerance  of  0.03  in.  or  0.2b  percent.  Con¬ 
stancy  of  water  depth  is  necessary  to  ensure  repeatability  of  test  conditions. 

Wave  Generation  and  Measurement 


62.  This  phase  is  controlled  by  one  process  called  TAPEM2.  Appendix  E 
outlines  the  procedures  and  inputs  required.  This  process  allows  the  reading 
of  a  control  signal  from  a  9T  magnetic  tape  containing  multiple  files  created 
in  the  control  signal  generation  phase  by  program  COMPONM  or  C0MP0NC4.  TAPEM2 
queries  the  user  for  the  length  of  time  to  run  a  test  and  automatically  turns 
oft  at  the  end  of  this  time. 

63.  The  length  of  time  data  is  collected  is  determined  by  the  user's 
response  to  the  query  for  measurement  delay  time  and  the  number  of  periods 
specified  in  Process  IDCAL  in  Header  2.  A  10-sec  time  period  was  allowed  to 
elapse  after  starting  of  tne  test  before  activating  the  DSWG.  This  allowed 
time  for  the  operator  to  turn  on  the  DSWG  if  no  assistance  was  available  and 
to  always  begin  at  the  same  time  tor  each  test.  For  a  2-sec  period,  6  cycles 
of  data  (12-sec  total  test  duration)  would  be  collected  after  the  DSWG  had 
been  running  tor  16  sec.  At  the  completion  of  data  collection,  the  DSWG  would 
continue  running  for  13  more  seconds  (i.e.  25  to  12  sec).  Based  on  results 
from  test  cases,  delay  times,  test  lengths,  and  collection  intervals  shown  in 
Table  9  for  the  five  wave  periods  were  used  in  all  tests.  Table  10  gives  ab¬ 
solute  times  from  start  to  finish  for  each  wave  period. 

64.  Prior  to  starting  a  test,  gages  were  zeroed  to  within  ±50  mv  using 
the  "Balance  R"  screw  adjustments  for  each  gage  on  the  instrument  console 
rack.  The  DSWG  is  activated  by  pressing  the  "Run"  program  button  on  the  MTS 
system  control  panel  of  the  wave  generator  control  console. 

65.  During  data  collection  it  is  important  to  prohibit  input/output 
(I/O)  from  other  users  as  this  may  interfere  with  accurate  generation  of  waves 
and  collection  of  data.  If  another  user  attempts  an  I/O  operation  (i.e. 
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logging  on,  editing,  file  manipulation,  compiling,  linking,  or  running  a  pro¬ 
gram,  etc.)  while  a  test  is  under  way,  the  control  signal  to  the  DSWG  is  mo¬ 
mentarily  interrupted  and  subsequently  accelerated  to  match  the  current  signal 
position.  The  data  collection  might  also  be  interrupted  because  the  system 
resources  are  momentarily  unavailable  for  data  collection  and  processing.  To 
prevent  this,  a  warning  message  is  sent  to  all  on-line  users  prior  to  taking 
data  to  alert  them  that  a  test  was  about  to  begin.  After  completing  a  test, 
the  "Stop"  program  button  on  the  MTS  system  control  panel  is  pressed  and  the 
"On  Line"  button  on  the  tape  drive  activated  for  the  next  test.  A  message 
sent  to  all  current  users  informs  them  that  testing  is  completed  and  they  may 
resume  normal  I/O  operations  until  further  notice. 

Wave  Analysis 

66.  After  the  data  have  been  collected,  analysis  takes  place  in  three 
stages.  Stage  1  is  a  preprocessing  stage  in  which  the  calibration 
coefficients  and  header  information  created  by  Process  IDCAL 

(DUAlrlDATA_STOjFOR018.DAT)  are  combined  with  data  collected  by  Process  TAPEM2 
(DLAO: [WAVEjACCEPT.DAT)  into  one  disk  file  in  the  directory  DUA1 : [DATA_ANAL] . 
Stage  2  is  the  zero-crossing  analysis  to  calculate  average  and  significant 
wave  period  and  height  for  model  and  prototype  (same  as  model  in  these  tests) 
waves.  A  harmonic  analysis  to  determine  the  percent  of  the  total  variance 
contained  in  a  fundamental  and  the  first  three  harmonics  is  calculated  in  the 
last  stage.  Command  procedure  files  for  each  of  these  three  stages  are  con¬ 
tained  in  Appendix  F.  Appendix  F  also  describes  the  operating  procedure  for 
each  stage. 
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PART  V:  TEST  RESULTS  AND  ANALYSIS 


67.  In  this  part,  results  from  measurements  of  wave  profiles,  periods, 
heights,  and  directions  for  both  linear  and  nonlinear  waves  are  presented  and 
discussed. 

Test  Case  Results 

68.  As  discussed  in  Part  III,  a  test  case  was  run  for  each  of  the  five 
wave  periods  to  verify  optimum  times  for  collecting  data.  Table  11  summarizes 
these  test  case  parameters  and  resulting  data  collection  intervals  selected. 

In  all  cases  a  stroke  of  1  in.  and  a  wave  direction  of  0  deg  were  used.  Tests 
were  run  for  a  long  duration  to  ascertain  the  degree  of  variation  in  wave  pro¬ 
file  with  time.  Travel  times  listed  are  theoretical  times  required  for  waves 
to  travel  from  the  DSWG  to  the  back  gage  row  (Timel)  and  to  reflect  off  the 
beach  and  return  to  the  back  gages  (Time2) .  (See  Part  ill  for  detailed  expla¬ 
nation  of  these  terms.)  Actual  data  collection  intervals  (repeated  here  for 
convenience  from  Tables  9  and  10)  were  determined  based  on  a  comparison  of 
these  values  with  wave  profiles  measured  in  the  back  gage  row  at  X  »  30  ft 
from  the  DSWG  (i.e.  gages  7,  8,  and  9).  For  gages  7,  8,  and  9  and  wave  peri¬ 
ods  of  0.75  and  1.50  sec.  Figures  I3a-c  and  14a-c,  respectively,  are  typical 
of  the  variations  in  wave  profiles  encountered.  Each  figure  has  a  two-line 
descriptive  title  at  the  top  for  documentation  purposes.  The  first  line  gives 
test  case  identification  code,  wave  direction  in  degrees,  DSWG  stroke  in 
inches,  and  number  of  cycles  of  the  wave  and  its  period  in  seconds.  The  sec¬ 
ond  header  line  gives  run  number  and  gage  number  for  the  test.  Note  the  in¬ 
crease  in  wave  profile  for  the  1.5U-sec  wave  period  after  approximately  60  sec 
elapsed  time.  This  type  of  variation  in  wave  profile  (averaged  for  all  three 
gages)  was  compared  with  the  theoretically  predicted  times  in  determining  op¬ 
timum  intervals  for  data  collection. 


Wave  Protile  Analysis 


69.  Representative  surface  elevation  time  series  for  each  of  the  five 
wave  periods  0.75,  1.00,  1.50,  2.U0,  and  3.00  sec  at  a  fixed  stroke  of  1  in. 
and  wave  direction  of  0  deg  are  shown  in  Figures  15a-e,  respectively.  The 
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Figure  15.  Wave  period  profile  (Sheet  1  of  3) 
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e.  Waves  of  3.00-sec  period 
Figure  15.  (Sheet  3  of  3) 
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measurement  location  waa  tha  cantar  of  tha  flrat  row  of  tha  measurement  area 
at  coordlnataa  X  -  10  ft,  Y  ■  45  ft  (l.a.  vava  gage  2).  Tha  amount  of  varia¬ 
tion  In  tha  wava  profile  with  tine  la  typical  of  that  observed  In  most  cases, 
with  aoaa  gagas  exhibiting  aore  and  some  gagaa  leas  variation.  Tha  measured 
wava  profiles  are  fairly  linear  for  these  low  wava  steepness  and  Goda's  non¬ 
linear  parameter  (NLP)  values.  Figures  16a-e  Illustrate  the  variation  of  wave 
profile  with  strokes  of  1,  3,  6,  9,  and  11  In.,  respectively,  for  a  fixed  pe¬ 
riod  of  3.00  sac  and  wave  direction  of  0  dag.  Tha  affect  of  Increasing  non¬ 
linearity  due  to  tha  binding  of  a  higher  harmonic  wava  component  for  strokes 
greater  than  3  In.  Is  evident  fron  these  figures. 

70.  Finally,  variation  In  wave  profile  with  the  five  wave  directions  0, 
15,  30,  45,  and  60  deg  is  Illustrated  In  Figures  17a-e,  respectively,  for  a 
fixed  wave  period  of  1.50  sec  and  a  stroke  of  3  in.  The  Increase  In  wave 
height  with  direction  Is  clearly  evident  In  these  figures,  between  15  and 
JO  deg,  the  Increase  Is  not  as  drastic  as  In  the  other  cases,  however. 

Wave  Period  Analysis 


Zero-crossing  analysis 

71.  Results  of  the  zero-crossing  wave  period  analysis  for  the  linear 
wave  cases  are  listed  In  Table  12  for  wave  directions  0,  15,  30,  45,  and 

60  deg.  The  aeasured  wave  period  versus  theoretical  value  is  given  for  each 
of  the  five  DSWG  stroke  wave  conditions.  Blank  spaces  In  the  tables  for  some 
wave  conditions  are  due  to  the  Inability  to  generate  a  particular  wave  condi¬ 
tion  because  of  wave  breaking.  Measured  average  values  are  the  average  for 
each  wave  period  over  all  strokes.  The  basin  response  factor  (BRF)  in  the 
last  column  Is  a  measure  of  the  efficiency  of  the  basin.  It  Is  the  ratio  of 
the  measured  average  to  the  theoretical  value  In  the  first  column.  The  BRF  Is 
also  used  In  the  paragraphs  which  follow  to  quantify  the  efficiency  of  the 
basin  relative  to  wave  height  and  wave  direction.  The  overall  BRF  for  all 
linear  wave  cases  at  0  deg  Is  99.65  percent.  Flguxe  18  shows  measured  versus 
theoretical  wave  periods  corresponding  to  Table  12  for  each  of  the  five  wave 
directions.  The  solid  curve  In  each  figure  Is  for  theoretical  values  and  the 
symbols  denote  different  wavemaker  strokes.  As  can  be  seen  from  the  plots  and 
the  BRF's,  agreement  Is  quite  good. 

72.  Table  13  summarizes  measured  average  wave  periods  from  Table  11  for 
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Figure  16.  Wave  profile  with  strokes  (Sheet  1  of  3) 
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Waves  of  60-deg  wave  direction 
Figure  17.  (Sheet  3  of  3) 
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Figure  18.  Measured  versus  theoretical  wave  periods, 
depth  -  1  ft  (Sheet  1  of  3) 
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Measured  Period.  Sec 
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e.  Waves  of  60  deg 
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Measured  Period,  Sec 


each  of  the  five  wave  directions.  An  overall  average  and  BRF  are  listed  in 
the  last  two  columns.  Figure  19  illustrates  excellent  agreement  between  mea¬ 
sured  and  theoretical  values  corresponding  to  an  average  BRF  of  99.6  percent 
for  all  strokes  and  directions. 


Theoretical  Period,  Sec 

Figure  19.  Measured  average  versus  theoretical  wave  period  for 
all  directions,  depth  -  I  ft 


Harmonic  analysis 

73.  As  discussed  previously  in  Part  III,  a  least-squares  harmonic  anal¬ 
ysis  was  performed  on  the  linear  wave  data  to  quantify  variance  contained  in 
the  fundamental  and  first  three  harmonics.  Table  14  lists  percent  variance  or 
total  energy  contained  in  the  first  harmonic  for  each  of  the  five  wave  direc¬ 
tions  0,  15,  30,  45,  and  60  deg.  In  addition  to  percent  variance,  Goda’s  NLP 
II  is  given  for  the  five  DSWG  stroke  and  wave  period  combinations  for  each 
wave  direction.  The  variation  of  total  energy  with  this  nondimensional  param¬ 
eter  is  shown  in  Figure  20  for  all  five  wave  directions.  Wave  conditions  with 
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Figure  20.  Harmonic  analysis  of  linear  waves 


n  values  less  than  0.2  have  approximately  85  percent  or  more  of  their  total 
energy  in  the  first  harmonic.  This  Indicates  a  reasonably  linear  or  sinusoi¬ 
dal  wave  profile. 


74.  Results  of  the  zero-crossing  wave  height  analysis  for  linear  wave 
cases  are  listed  in  Table  15  for  the  five  wave  directions  0,  15,  30,  45,  and 
60  deg.  In  the  table  measured  wave  height  is  compared  with  the  theoretical 
value  (see  Part  III)  for  each  wave  period  and  stroke  combination.  The  value 
of  maximum  possible  wave  height  (i.e.  prebreaking  wave  height)  for  each  wave 
period  is  repeated  for  reference  in  the  second  column  of  each  table.  Maximum 
measured  values  are  discussed  in  paragraph  77.  Again,  blank  spaces  in  the 
measured  columns  are  because  it  is  not  possible  to  generate  waves  with  these 
conditions  without  breaking  occurring.  As  explained  In  Part  III,  the  effect 
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of  directionality  on  wave  height  is  to  increase  its  value  relative  to  the  two- 
dimensional  height  of  a  wave  traveling  perpendicular  to  the  DSWG  (i.e.  0-deg 
wave  direction).  Thus,  as  the  wave  direction  increases,  the  number  of  waves 
which  can  be  generated  tend  to  decrease  in  "stairstep"  fashion  as  the  wave 
period  and  stroke  increase. 

75.  Plots  of  the  measured  versus  theoretical  wave  heights  for  each  of 
the  five  strokes  tested  are  shown  in  Figure  21  for  wave  directions  corre¬ 
sponding  to  Table  14.  The  solid  curves  in  each  figure  represent  theoretical 
values  and  the  symbols  denote  different  wavemaker  strokes.  The  top  curve  in 
each  section  of  the  figure  is  the  "breaking  wave  limit  curve,"  based  on  Equa¬ 
tion  5.  The  theoretical  curves  were  plotted  from  calculations  based  on  only 
the  five  discrete  wave  periods  tested.  Thus,  they  are  not  as  smooth  as  they 
would  have  been  had  they  been  plotted  using  a  continuous  wave  period  array. 

76.  From  the  table  and  plots,  it  is  evident  that  linear  wave  theory 
usually  overpredicts  wave  height  for  the  DSWG  basin.  The  only  exception  ap¬ 
pears  to  be  the  3-sec  period,  11-in.  stroke  combination  for  wave  directions 
less  than  30  deg.  The  average  BRF  (see  paragraph  71)  for  all  period  and 
stroke  combinations  at  0  deg  is  87.7  percent.  The  average  BRF  over  all  wave 
cases  is  79.2  percent.  The  general  trend  is  for  the  BRF  to  Increase  with  in¬ 
creasing  stroke  for  a  fixed  wave  direction  and  to  decrease  with  wave  direction 
for  all  stroke  combinations. 

Wave  breaking  tests 

77.  To  verify  maximum  possible  wave  heights,  a  series  of  three  waves 
was  run  tor  each  wave  period  to  bracket  the  breaking  phenomenon.  All  tests 
were  run  at  a  wave  direction  of  0  deg.  Table  16  and  Figure  22  Illustrate  the 
results  from  this  series  of  tests.  The  results  are  in  agreement  with  those 
reported  previously;  namely,  the  wave  heights  are  slightly  overpredicted  for 
the  DSWG  basin.  Measured  wave  heights  near  breaking  always  occurred  below 
predicted  values,  especially  for  wave  periods  greater  than  1.5  sec.  Required 
strokes  always  were  larger  than  predicted  theoretical  values.  For  the  3-sec 
wave  period,  it  was  not  possible  to  determine  the  rccuracy  of  linear  theory 
because  a  stroke  larger  than  the  maximum  DSWG  stroke  of  12  in.  is  required. 
Constancy  of  wave  height 

78.  To  evaluate  constancy  of  wave  height  within  the  measurement  area, 
the  measured  wave  height  at  each  of  the  nine  gages  was  examined  for  linear 
control  signals  with  strokes  of  3  in.,  periods  of  1.00  and  1.50  sec,  and  wave 
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Figure  21.  Measured  versus  theoretical  wave  heights, 
depth  ■  1  ft  (Sheet  1  of  3) 
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Figure  21.  (Sheet  2  of  3) 
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Figure  21.  (Sheet  3  of  3) 
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Figure  22.  Wave  breaking  tests 

directions  of  0,  15,  30,  and  45  deg.  Table  17  lists  the  theoretical  (based  on 
the  three-dimensional  height-to-stroke  H/S  ratio),  measured,  average  H  , 
and  mean  variation  wave  height  AH  for  each  of  the  eight  cases  for  both  small 
and  Jordan  wave  gages.  The  purpose  of  the  gage  comparisons  is  to  quantify 
what  effect,  if  any,  the  higher  resolution  of  the  small  gages  has  on  the  mea¬ 
sured  values.  Corresponding  measured  wave  heights  with  a  period  of  1.00  sec 
for  both  small  and  Jordan  wave  gages  are  plotted  in  Figure  23  for  wave  direc¬ 
tions  0,  15,  30,  and  45  deg.  Similarly,  the  cases  with  a  period  of  1.5  sec 
are  shown  in  Figure  24  for  their  wave  directions.  The  data  tor  the  small  wave 
gages  are  shown  with  a  dashed  line  and  that  ot  the  larger  range  Jordan  wave 
gages  with  a  solid  curve. 

79.  Although  some  measured  wave  heights  appear  to  vary  significantly 
among  gages,  AH  calculated  for  the  array  of  nine  gages  was  always  within 
tolerance  limits  suggested  by  Sand  (1979).  (See  Part  III.)  These  values  fol¬ 
lowed  the  same  pattern  he  observed  of  increasing  with  increasing  wave  direc¬ 
tion.  The  one  exception  was  for  30  deg,  where  AH  seemed  to  dip  slightly  in 
all  cases  except  for  the  1.00-sec  Jordan  gages.  The  1,50-sec  wave  periods 
seemed  to  give  less  wave  height  variation  than  corresponding  1.00-sec  wave 
periods,  except  for  a  wave  direction  of  45  deg.  For  the  l. 00-sec  waves,  there 
there  is  a  large  scatter  of  measured  points  about  the  center  line  of  the  mea¬ 
surement  area  for  different  directions. 
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Figure  23.  Constancy  of  wave  heights  (S10  means 
small  gage,  X  -  10  ft;  J10  means  Jordan  gage, 

X  -  10  ft,  etc.)  (Continued) 
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Distance  Parallel  to  DSWC,  Ft 


a.  Waves  of  1.50  sec  and  0  deg 
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b.  Waves  of  1.50  sec  and  15  deg 
Figure  24.  Constancy  of  wave  heights  (Continued) 


80.  To  ascertain  repeatability,  the  tests  were  run  twice  for  small 
gages  at  a  wave  direction  of  0  deg.  Measured  wave  heights  for  each  gage  were 
very  consistent,  differing  no  more  than  0.5  percent  in  AH  between  runs. 

81.  The  major  cause  of  variation  is  probably  higher  harmonics  (espe¬ 
cially  free  second)  and  wave  reflections  within  the  asymmetric  basin.  Wave 
absorption  characteristics  of  the  basin,  as  described  in  Part  II,  are  quite 
reasonable.  The  combination  of  all  sources  of  variability,  as  described  in 
Part  III,  contributes  to  overall  measured  wave  height  variability. 

82.  An  additional  explanation  of  some  observed  variation  is  the  toler¬ 
ance  on  wave  gage  calibrations.  The  maximum  deviation  of  each  wave  gage  at 
any  one  step  in  the  11-step  calibration  process  is  a  good  indicator  of  this 
tolerance  (see  Appendix  Tables  D3  and  D4) .  The  average  values  of  maximum  de¬ 
viation  for  small  and  Jordan  gages  are  0.015  in.  (0.47  percent  of  full  scale 
of  3.25  in.)  and  0.145  in.  (1.45  percent  of  full  scale  of  10  in.),  respec¬ 
tively.  The  Jordan  gages'  maximum  deviation  values  are  approximately  10  times 
larger  than  those  for  the  small  gages.  Relative  to  full  scale  of  calibration, 
the  difference  between  the  two  is  only  about  three  times  for  the  Jordan  gage. 

Wave  Direction  Analysis 


Integer  paddle  assumption 

83.  As  explained  in  Part  III,  an  integer  number  of  paddles  was  used  to 
determine  wave  directions  in  the  control  signal.  Thus,  except  for  the  wave 
direction  of  0  deg,  measured  wave  directions  do  not  match  desired  wave  direc¬ 
tions  of  15,  30,  45,  and  60  deg  exactly  but  are  reasonably  close.  Measured 
wave  directions  are  compared  with  input  wave  directions  rather  than  with 
desired  wave  direction.  The  amount  of  variation  increases  for  directions 
greater  than  30  deg  and  wave  periods  less  than  1  sec.  The  average  BRF  (see 
paragraph  71)  for  the  wave  direction  of  0  deg  is  99.9  percent.  The  overall 
BRF  for  all  directions,  strokes,  and  period  combinations  was  98.81  percent, 
indicating  a  high  probability  of  generating  a  wave  with  the  desired  wave 
direction. 

84.  Table  18  lists  measured  versus  theoretical  wave  directions  fox  the 
five  desired  wave  directions  0,  15,  30,  45,  and  60  deg.  The  format  of  this 
table  is  similar  to  that  for  periods  and  wave  heights  presented  earlier. 
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Figure  25  (corresponding  to  Table  18)  illustrates  these  relationships  versus 
wave  period  for  each  desired  wave  direction. 

85.  Table  19  lists  measured  average  values  (eighth  column  from 
Table  18)  for  each  desired  wave  direction  versus  wave  period.  The  maximum 
wave  direction  obtainable  with  an  Integer  number  of  paddles  for  each  wave  pe¬ 
riod  is  repeated  in  the  second  and  third  column  (see  Part  III).  Finally,  the 
excellent  agreement  between  measured  and  theoretical  wave  directions  is 
plotted  in  Figure  26.  The  solid  lines  are  theoretical  values,  and  symbols 
denote  the  average  for  each  direction  over  all  period  and  stoke  combinations. 
Again,  the  top  curve  is  the  breaking  wave  limit  curve  from  Equation  5.  The 
worst  fit  is  for  the  3-sec  wave  period  at  60  deg  where  the  average  BFR  is 
95.1  percent. 

Exact  paddle  assumption 

86.  In  order  to  verify  the  ability  of  the  DSWG  to  generate  wave  direc¬ 
tions  requiring  a  noninteger  number  of  paddles,  a  series  of  three  tests  was 
run  with  a  fixed  wave  period  of  1.5  sec  and  a  desired  wave  angle  of  approxi¬ 
mately  60  deg.  The  second  test  was  identical  to  the  case  with  an  Integer  num¬ 
ber  of  paddles.  The  other  two  cases  were  designed  to  bracket  the  desired 
value  (i.e,  60  deg)  with  a  noninteger  number  of  paddles  to  generate  waves  with 
directions  above  and  below  60  deg.  Based  on  the  values  measured  and  the  BRF's 
calculated  (Table  20) ,  wave  directions  near  ±90  deg  can  be  generated  at  any 
wave  period,  provided  more  than  two  paddles  are  used.  For  wave  periods  less 
than  1  sec,  however,  the  ability  to  generate  waves  near  ±90  deg  is  somewhat 
limited  since  less  than  three  paddles  are  required.  As  pointed  out  in 

Part  ill,  a  more  reasonable  upper  limit  on  0.75-sec  waves  is  70  deg. 

Nonlinearity  Effects 

87.  To  ascertain  the  importance  of  nonlinear,  second-order  control  the¬ 
ory,  a  series  of  16  cnoldal  waves  was  created  and  tested  in  the  DSWG  basin. 
These  waves  consisted  of  combinations  of  three  periods  at  1.5,  2.0,  and 

3.0  sec  and  three  generator  strokes  of  6,  9,  and  11  in.  for  each  of  two  wave 
directions  of  0  and  30  deg.  It  is  not  possible  to  generate  a  cnoldal  wave 
with  a  3-sec  period  and  11-in.  stroke  in  a  water  depth  of  1  ft.  The  control 
and  feedback  signals  were  monitored  to  verify  accuracy  of  the  control  signal 
in  form,  period,  and  amplitude. 
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Wav«  Direction,  0«9  Wov«  Dlr«etton, 


a.  Waves  of  0  deg 


b.  Waves  of  15  deg 

Figure  25.  Measured  versus  theoretical  wave  directions, 
depth  ■  1  ft  (Sheet  1  of  3) 


c.  Waves  of  30  deg 
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Figure  25.  (Sheet  2  of  3) 
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Wove  Direction. 


e.  Waves  of  60  deg 


Figure  25.  (Sheet  3  of  3) 
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Figure  26.  Measured  average  versus  theoretical  wave 
directions  for  all  directions,  depth  ■  1  ft 


wave  profile  analysis 

»8.  Effect  of  wave  period  on  linear-  and  nonllnear-generated  wave  pro¬ 
files  Is  Illustrated  In  Figure  27.  Figures  27a-c  show  wave  profiles  for  lin¬ 
ear  waves  measured  at  the  center  of  the  first  row  of  the  measurement  area 
(l.e.  X  -  10  ft,  Y  *  45  ft)  for  a  b-in.  stroke  and  0-deg  wave  direction  for 
the  three  wave  periods  1.5,  2.0,  and  3.0  sec,  respectively.  Figures  27d-f 
illustrate  comparable  wave  profiles  for  nonlinear  waves.  The  linear  control 
signal  profiles  Illustrate  the  binding  of  higher  harmonic  waves  to  form  a  non¬ 
linear  wave.  Nonlinear  control  signal  profiles  for  the  3.00-sec  period  (Fig¬ 
ure  27f)  Illustrate  a  possible  mismatch  In  the  control  signal  as  an  FHW  is 
evident  along  the  front  side  of  the  crests. 

#9.  Figure  28  illustrates  the  effect  of  wavemaker  stroke  on  linear-  and 
nonllnear-generated  wave  profiles.  Figures  28a-c  show  measured  results  for 
linearly  generated  control  signals  for  a  fixed  3.00-sec  period  and  wave  direc¬ 
tion  of  0  deg  for  the  three  wavemaker  strokes  of  6,  9,  and  II  in.,  respec¬ 
tively.  Again,  Figures  28d-e  show  analogous  nonlinear  wave  profiles.  The 
linear  control  signal  profiles  are  again  nonlinear  In  shape.  The  nonlinear 
control  signal  profiles  Indicate  a  strong  FHW  for  both  6-  and  9-in.  wavemaker 
strokes. 
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b.  Linear  waves  of  2.00-sec  period 
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Figure  27.  Wave  profile,  effect  of  wave  period 
(Sheet  1  of  3) 
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Figure  27.  (Sheet  2  of  3) 
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f.  Nonlinear  waves  of  3.00-sec  period 
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Linear  waves  of  6-in.  stroke 
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Linear  waves  of  9-in.  stroke 


Figure  28.  Wave  profile,  effect  of  waveraaker  stroke 
(Sheet  1  of  3) 
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d.  Nonlinear  waves  of  6-in.  stroke 
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Figure  28.  (Sheet  3  of  3) 


90.  Finally,  Figure  29  shows  the  effect  of  wave  direction  on  wave  pro¬ 
files  for  linear  and  nonlinear  control  signals.  Figures  29a-b  illustrate  this 
directional  eftect  for  linear  wave  cases  for  a  fixed  period  of  1.50  sec, 
stroke  of  6  in.,  and  wave  directions  of  0  and  30  deg,  respectively.  Fig¬ 
ures  29c-d  show  the  corresponding  effect  on  nonlinear  wave  signals.  The 
three-dimensional  effect  of  directionality  is  evident  in  these  figures. 
Transformation  of  the  linearly  generated  wave  to  a  nonlinear  shape  Indicates 
the  need  to  use  nonlinear  control  signals  for  these  particular  combinations  of 
wave  parameters. 

Wave  period  analysis 

91.  A  comparison  of  the  measured  and  theoretical  nonlinear  wave  periods 
for  0-  and  30-deg  wave  directions  is  given  in  Table  21.  The  overall  agreement 
is  quite  good,  with  BRF's  of  99.7  percent  for  both  the  0-  and  30-deg  wave 
directions . 

92.  Goda  (1983)  verified  the  applicability  of  his  nonlinear  parameter 
II  for  describing  the  nonlinearity  of  water  waves  from  deep  to  shallow  water 
using  regular  laboratory  waves  and  strokes  third-order  and  second-order  cnoi- 
dal  wave  theory.  Nonlinear  cnoidal  waves  can  be  decomposed  into  the  relative 
magnitudes  of  their  Fourier  components  of  wave  profiles  using  harmonic  analy¬ 
sis.  The  more  nonlinear  the  wave  profile,  the  larger  the  percent  variance  in 
the  higher  harmonics  relative  to  the  first  harmonic.  The  average  percent 
variance  in  the  first  harmonic  for  nine  gages  versus  Goda's  nonlinear  parame¬ 
ter  II  for  linear  and  nonlinear  waves  is  presented  in  Table  22  and  Figure  30. 
For  both  linear  and  nonlinear  waves,  this  value  decreases  with  increasing  II  . 
There  appears  to  be  a  crossover  at  a  certain  value  of  II  where  the  nonlinear 
control  signal  better  represents  the  desired  waveform  (i.e.,  the  variance  in 
the  first  harmonic  is  larger  for  the  nonlinear  wave  relative  to  the  linear 
wave).  Based  upon  the  limited  data  available,  this  appears  to  be  approxi¬ 
mately  a  II  of  0.27  to  0.29. 

Wave  height  analysis 

93.  Table  23  shows  measured  versus  theoretical  wave  heights  for  non¬ 
linear  waves  for  wave  directions  0  and  30  deg  (assuming  that  a  wavemaker 
stroke  of  6,  9,  or  11  in.  will  generate  a  wave  height  of  comparable  value). 

In  all  cases,  the  measured  wave  heights  were  much  smaller  than  the  assumed 
values.  Thus,  the  BRF  factors  reported  here  are  representative  of  more  than 
just  the  DSWG  basin  efficiency  as  defined  for  linear  waves,  but  also  include 
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Figure  29.  Wave  profile,  effect  of  wave 
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Goda's  Nonlinear  Parameter 

Figure  30.  Linear  versus  nonlinear  harmonic 
analysis  results 


the  effect  of  the  wave  generator  transfer  function.  The  measured  wave  heights 
decrease  with  increasing  period  for  a  fixed  value  of  DSWG  stroke,  just  as  for 
linear  waves.  It  also  decreases  with  increasing  stroke  for  a  fixed  wave  pe¬ 
riod.  Breaking  occurred  in  some  of  the  cases  reported.  This  resulted  in  a 
smaller  wave  height  than  would  be  present  had  breaking  and  subsequent  loss  ot 
energy  not  occurred.  The  purpose  of  reporting  these  values  is  to  give  some 
indication  of  the  wave  height  magnitude  and  where  wave  breaking  becomes  a 
serious  consequence  for  nonlinear  waves. 

Wave  direction  analysis 

94.  Table  24  lists  measured  versus  theoretical  wave  directions  for  de¬ 
sired  wave  directions  of  0  and  30  deg.  Agreement  is  excellent  tor  the  0-deg 
cases.  Measured  values  were  considerably  higher  for  a  desired  wave  angle  of 
30  deg  for  nonlinear  waves  relative  to  the  linear  waves  (see  Table  18).  The 
BRF  for  a  wave  direction  of  0  deg  is  99.9  percent  and  82.9  percent  for  30  deg. 
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PART  VI:  SUMMARY  AND  RECOMMENDATIONS 


95.  The  directional  spectral  wave  generator  (DSWG)  is  a  unique  re¬ 
source  for  the  study  of  natural  sea  states  in  a  laboratory  environment.  Its 
unique  features  Include  size,  modular  design,  portability,  method  of  paddle 
connection  and  displacement,  and  electric  motor  power.  More  economical  and 
efficient  design  of  coastal  structures  will  result  because  of  the  use  of  the 
DSWG  facility  in  site-specific  and  research  studies. 

96.  A  series  of  111  linear  and  nonlinear  control  signals  were  gener¬ 
ated  to  quantify  performance  characteristics  of  the  DSWG  basin.  Measurements 
were  made  at  nine  locations  within  the  basin  by  resistance-wire  wave  gages. 
Comparisons  of  measured  versus  predicted  wave  profiles,  periods,  heights,  and 
directions  Indicate  the  range  of  applicability  ot  linear  control  signal 
theory. 

97.  The  linear  control  signals  consisted  of  combinations  of  five  wave 
periods  (i.e.  0.75,  1.00,  1.50,  2.00,  and  3.00  sec),  five  wavemaker  strokes 
(i.e.  1,  3,  6,  9,  and  11  in.),  and  five  wave  directions  (i.e.  0,  15,  30,  45, 
and  60  deg)  representative  of  model  test  conditions  for  monochromatic  waves. 

Of  the  125  possible  combinations  (i.e.  five  periods  by  five  strokes  by  five 
directions) ,  only  80  were  actually  tested  because  certain  combinations  of 
period,  stroke,  and  direction  are  impossible  to  generate  without  wave  breaking 
occurring. 

98.  A  basin  response  factor  (BRF)  was  calculated  to  indicate  the  abil¬ 
ity  of  the  basin  to  accurately  reproduce  theoretically  predicted  wave  parame¬ 
ters.  It  is  the  ratio  of  the  measured  average  value  for  all  nine  gages  to  the 
theoretical  value.  The  measured  wave  periods  show  excellent  agreement  with 
the  predicted  values.  The  average  BRF  for  all  wave  strokes  and  directions  is 
99.6  percent. 

99.  A  least-squares  harmonic  analysis  was  performed  to  quantify  the 
variance  (total  energy)  contained  in  the  fundamental  harmonic  component.  The 
larger  this  value,  the  more  linear  the  actual  waveform.  Goda's  nonlinear 
parameter  II  also  was  calculated.  Wave  conditions  with  II  values  less  than 
0.2  have  approximately  85  percent  or  more  of  their  total  energy  (variance)  in 
the  first  harmonic.  This  Indicates  a  reasonably  linear  wave  profile. 

100.  The  height-to-stroke  ratios  predicted  by  linear  wave  theory  gener¬ 
ally  overpredict  the  measured  wave  heights  for  different  wave  period,  stroke. 
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and  direction  combinations.  The  only  exception  was  for  3-sec  period,  11-in. 
stroke  waves  of  direction  less  than  30  deg.  The  average  BRF  for  all  period 
and  stroke  combinations  at  0  deg  is  87.7  percent,  decreasing  to  79.2  percent 
for  all  directions.  The  general  trend  is  for  the  BRF  to  increase  with  in¬ 
creasing  stroke  for  a  fixed  wave  direction  and  to  decrease  with  wave  direction 
for  all  stroke  combinations.  A  major  cause  for  lower  than  predicted  wave 
heights  is  the  fact  that  the  DSWG  is  not  sealed  at  the  ends  and  along  the  bot¬ 
tom.  This  results  in  a  loss  of  energy  from  the  wave  field  as  waves  are  gener¬ 
ated  and  a  corresponding  loss  in  generated  wave  height.  Measured  wave  heights 
for  different  wave  directions  did  not  always  agree  well  with  predicted  values 
using  first-order  theory.  Additional  research  on  second-order  effects  of 
directionality  on  height-to-stroke  ratios  would  be  beneficial. 

101.  Constancy  of  mean  wave  height  was  calculated  for  several  cases 
with  wave  periods  of  1.00  and  1.50  sec.  Although  some  measured  wave  heights 
appear  to  vary  more  than  expected,  the  variation  in  mean  wave  height  AH  was 
always  less  than  the  15  to  25  percent  tolerance  limits  suggested  by  Sand 
(1979).  The  value  of  AH  always  followed  the  same  pattern  observed  by  Sand 
of  increasing  with  increasing  wave  direction.  The  most  probable  sources  of 
this  variation  are  the  interaction  of  free  and  bound  higher  harmonics,  wave 
reflections,  end  diffraction,  Benjamin-Feir  type  of  effect  due  to  low  fre¬ 
quency  standing  cross  waves,  asymmetry  of  the  basin,  variable  bottom  gap,  and 
wave  gage  calibration  tolerances.  Additional  research  is  needed  in  these 
areas,  especially  prediction  and  control  of  higher  harmonics,  improved  wave 
absorption  for  the  beaches  and  side  walls,  and  effect  of  calibration  toler¬ 
ances  on  repeatability  of  wave  gages. 

102.  The  design  of  connections  between  paddles  of  the  DSWG  allows  gen¬ 
eration  of  waves  of  a  directional  range  approaching  ±90  deg  for  all  wave  peri¬ 
ods  larger  than  1.00  sec.  For  periods  less  than  this,  the  maximum  wave  direc¬ 
tion  is  reduced  due  to  the  requirement  for  at  least  two  paddles  per  wavelength 
along  the  wave  generator.  A  more  reasonable  upper  limit  for  0.75-sec  waves 
using  slightly  more  than  two  paddles  is  about  70  deg.  The  agreement  between 
predicted  and  measured  wave  directions  was  excellent  for  all  linear  wave 
cases.  The  average  BRF  for  a  wave  direction  of  0  deg  is  99.9  percent,  de¬ 
creasing  to  98.8  percent  for  all  cases.  Thus,  the  snake  principle  works  very 
well  in  predicting  wave  directions  in  the  DSWG  basin.  Since  wave  directions 
were  calculated  manually  for  this  study,  it  is  recommended  that  a 
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cross-correlation  technique  be  Implemented  for  future  projects.  An  option 
should  be  included  to  calculate  the  phase  lag  between  pairs  of  wave  gages  and 
to  use  a  specific  gage  pair  or  the  average  of  all  possible  pairs. 

103.  Certain  wave  period  and  height  combinations  produce  waves  which 
tend  to  deviate  from  the  linear  wave  profile.  These  waves  become  more  non¬ 
linear  with  higher,  sharper  crests  and  broader,  flatter  troughs  due  to  a 
second-order  effect  of  a  bounded  (or  locked)  higher  harmonic  wave  component 
(BHW) .  If  linear  control  signal  theory  is  used  for  these  waves,  boundary  con¬ 
ditions  at  the  wavemaker  cannot  be  satisfied  for  the  second-order  BHW  (i.e., 

a  mismatch  in  wavemaker  motion  with  the  required  water  particle  velocities) 
resulting  in  the  creation  of  spurious  or  free  higher  harmonic  component 
waves  (FHW) .  These  FHW's  are  a  primary  cause  of  the  variations  in  wave 
height,  discussed  earlier,  due  to  alternate  cancellation  and  reinforcement 
with  the  BHW.  Thus,  to  ascertain  the  Importance  of  nonlinear,  second-order 
control  theory,  a  series  of  16  cnoidal  waves  was  created  and  tested  in  the 
DSWG  basin. 

104.  Comparisons  of  several  linear  wave  profiles  revealed  the  binding 
of  a  BHW  to  form  a  nonlinear  waveform.  The  nonlinear  profiles  for  3.00  sec 
illustrate  a  possible  mismatch  of  the  control  signal  as  an  FHW  is  evident 
along  the  front  side  of  the  crests. 

105.  The  wave  periods  for  nonlinear  waves,  as  in  the  linear  case, 
showed  good  agreement  with  predicted  values.  The  average  BRF's  for  both  the 
0-  and  30-deg  cases  were  99.7  percent. 

106.  A  comparison  of  harmonic  analysis  results  between  several  linear 
and  nonlinear  cases  indicates  that  nonlinear  control  signals  might  better  rep¬ 
resent  desired  waveforms  around  a  11  value  of  0.27  to  0.29.  Further  research 
into  this  relationship  with  Gods' s  nonlinear  parameter  is  recommended. 

107.  For  nonlinear  waves,  the  height-to-stroke  ratios  were  not  consid¬ 
ered  in  the  initial  determination  of  required  stroke  (as  was  done  for  linear 
wave  cases).  Thus,  measured  wave  heights  were  quite  a  bit  smaller  than  as¬ 
sumed  values,  and  the  BRF's  comprise  height-to-stroke  ratios  in  addition  to 
basin  efficiency.  Therefore,  it  is  not  valid  to  compare  them  with  comparable 
linear  BRF's.  Trends  are  the  same,  however,  as  for  linear  waves. 

108.  The  measured  versus  theoretical  wave  directions  for  directions  of 
0  deg  agreed  very  well  with  a  BRF  of  99.9  percent.  For  the  30-deg  wave  direc¬ 
tion,  however,  the  nonlinear  waves  do  not  show  as  good  agreement  as  the  linear 
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waves.  The  BRF  was  82.9  percent  with  the  measured  values  3  to  5  deg  above  the 
corresponding  linear  values. 

109.  Tests  with  nonlinear  control  signals  were  not  intended  to  be  ex- 
haustlve.  Additional  tests  should  be  run  to  better  quantify  applicable  ranges 
of  use  relative  to  linear  control  theory.  Also,  further  study  of  the  theoret¬ 
ical  development  and  verification  of  the  height-to-stroke  ratio  should  be 
pursued . 

110.  The  DSWG  basin  bathymetry  is  relatively  constant  over  the  90-  by 
70-ft  area  in  front  of  the  DSWG.  The  maximum  total  variation  is  1.08  in.  The 
highest  point  is  0.51  in.  (4.3  percent  In  1  ft  of  water)  and  the  lowest  point 
is  -0.57  in.  (4.8  percent  in  1  ft  of  water),  and  they  are  located  in  opposite 
corners  of  the  basin.  The  standard  deviation  for  the  basin  is  0.18  in.  For 
the  20-  by  20-ft  measurement  area,  the  maximum  total  variation  is  0.35  in. 

The  high  point  is  0.08  in.  (0.7  percent  in  1  ft  of  water)  and  the  low  point  is 
-0.27  in.  (2.3  percent  in  1  ft  of  water). 

111.  As  a  consequence  of  the  finite  width  of  wavemaker  paddles,  small 
undesired  disturbances  known  as  spurious  waves  are  generated  in  addition  to 
the  main  monochromatic  wave.  Since  the  DSWG  is  driven  at  the  joints  between 
paddles  rather  than  at  individual  paddles,  the  amplitude  of  any  spurious  waves 
generated  is  greatly  attenuated.  Spurious  wave  generation  was  not  a  serious 
problem  in  this  study  due  to  the  wave  generator  design  and  period  range 
tested. 
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Table  1 


ADACS  Computer  Specifications 

Specification 

Description 

Computer 

DEC  VAX  11/750 

Location 

Wave  Processes  Branch,  WES 
Bldg.  6006 

Operating  system 

VAX/VMS  Version  4.2 

Operating  mode 

Interactive,  real  time, 
prioritized  multitasking 

Language 

Fortran  77 

Memory 

4  MB 

Disk  storage 

577  MB  fixed 

10  MB  removable 

A/D  channels 

128  multiplexed, 
single-ended 

D/A  channels 

61  Preston,  IEEE  488 

4  DEC 

Tape  drives 

2  @  125  IPS,  800/1600  BPI 

Plotter 

Versatec  V-80 

Terminals 

VT  100  compatible 

Tektronix  4014,  4114-A 


Printer 


600  LPM 


138  columns 


Table  2 

Maximum  Wave  Direction  Versus  Wave  Period 


Exact  Values _  _ Integer  Values 


Period 

sec 

Wavelength 

ft 

Number 

Paddles 

Offset 

deg 

Direction 

deg 

Number 

Paddles 

Offset 

deg 

Direction 

deg 

0.75 

2.82 

1.88* 

191.5 

90 

3 

120 

38.8 

1.00 

4.52 

3.01 

119.6 

90 

4 

90 

48.9 

1.50 

7.73 

5.15 

69.9 

90 

6 

60 

59.2 

2.00 

10.77 

7. 18 

50.1 

90 

8 

45 

63.8 

3.00 

16.63 

11.09 

32.5 

90 

12 

30 

67.5 

*  Impossible  to  generate,  need  a  minimum  of  two  paddles. 


Table  3 

Number  of  Spurious  Waves  Based  on  Main  Wave  Parameters 


Period 

Wavelength 

Direction 

Number  of 

sec 

ft 

deg 

Waves 

0.25 

0.32 

0-5 

7+ 

0.50 

1.28 

0-8 

2 

8-25 

1 

0.75 

2.82 

0-60 

0 

60  -  68 

1 

1.00 

4.52 

0-90 

0 

Note:  Paddle  width 

B  -  1.5  ft;  water  depth 

h  -  1.0  ft. 

Table  4 

Wave  Properties  Examined 


_ Function _ 

Waveform 

Transfer 

Spurious  waves 
Wave  direction 
Waveform  nonlinearity 
Measurement  times 


_ Wave  Property _ 

Wavelength 
Celerity 
Group  speed 
Dimensionless  depth 

Height-to-stroke  ratio 
Two-dimensional  wave  height 
Three-dimensional  wave  height 
Maximum  prebreaking  wave  height 

Number  of  waves 
Minimum  wave  period 

Paddle  offsets 
Predicted  wave  direction 

Steepness  ratio 

Goda's  nonlinearity  parameter 

Beginning  delay  time 
Ending  delay  time 


AD-A182  142 
UNCLASSIFIED 


io 


SEARCH 
'f/C  28/4 


2/2 


NL 


Table  5 

Description  of  Subroutines  in  Program  MONOSUMMARY 


Name 

INPUT 

CASE 

HUNT 

SPEED 

KHCALC 

GROUP 

ANGLE2 

SPURIOUS 

HSRATIO 

HEIGHT 

HTHETA 

HBREAK 
STEEP 
NLP ARAM 
DELAY 

REFLECT 

HD120 


_ _ Description _ 

Queries  user  for  input  parameters 

Case  or  run  number  for  Individual  test  cases 

Hunt's  method  for  wavelength,  see  CETN-I-17* 

Vave  celerity  or  speed 
Nondimens ional  water  depth  ratio 
Group  speed 

Offset  angle  for  integer  number  of  paddles  to  produce  desired 
wave  direction  angle 

Minimum  wave  period  below  which  spurious  waves  will  be  generated 

Two-dimensional  wave-height-to-stroke  transfer  function 

Two-dimensional  wave  height,  same  as  three-dimensional  for  zero 
wave  direction 

Three-dimensional  wave  height  based  on  effects  of  wave  direction 
angle 

Maximum  breaking  wave  height  for  laboratory  waves 

Wave  steepness  ratio,  measure  of  nonlinearity  of  waveform 

Goda's  nonlinear  parameter 

Delay  time  prior  to  data  collection,  time  for  specified  number  of 
waves  to  reach  back  gage  row 

Trip  time  for  specified  wave  to  travel  from  wave  generation  to 
back  wall  and  return  to  back  gage  row 

Maximum  breaking  wave  height  based  on  120-deg  maximum  crest  angle 


*  US  Army  Engineer  Waterways  Experiment  Station,  Coastal  Engineering  Research 
Center  (1985). 


Table  6 

Input  Variables  for  Program  MONOSUMMARY 


NTH ETA  Number  of  wave  direction  angles  Maximum  number  ■  25 

THETA  Array  of  wave  angles,  deg  Enter  NTHETA  values 

NSTROKE  Number  of  generator  strokes  Maximum  number  ■  25 

STROKE  Array  of  generator  strokes.  In.  Enter  NSTROKE  values 

NPERIOD  Number  of  wave  periods  Maximum  number  *  25 

PERIOD  Array  of  wave  periods,  sec  Enter  NPERIOD  values 

H  Water  depth,  ft 

DISTNC1  Distance  to  back  gage  row,  ft  Used  in  delay  time  calculation 

NWAVE  Number  of  waves  to  pass  before  Used  in  delay  time  calculation 

testing 

DISTNC2  Trip  distance  to  back  wall  and  Maximum  time  before  reflected 

return  to  back  gage  row  energy  could  return 

TITLE2  Descriptive  title  for  project  .LE.  40  characters 

TITLE3  Descriptive  title  for  line  3  .LE.  40  characters 

SERIES  Test  series  description  1  character 

(i.e.  M  •  Monochromatic) 

MONTH  Month  of  test  1  character 

(i.e.  2  -  Feb) 


LOCATION  Test  location  1  character 

(i.e.  0  *  Initial  setup) 


Table  7 

Output  Variables  for  Program  MONO SUMMARY 


Name  _ Description _  _ Comment 


HDMAX 

Maximum  breaking  height 

Assumes  H/D  “  0.78 

NAME 

Case  or  run  number  Identification 

6  characters 

T 

Period,  sec 

L 

Shallow-water  wavelength,  ft 

Assumes  linear  theory 

C 

Wave  speed  or  celerity,  fps 

CG 

Group  speed,  fps 

OFFSET 

Offset  phase  angle  corresponding 
to  desired  wave  angle,  deg 

See  Part  III,  para  32 

DIR 

Wave  angle  generated,  deg 

See  Part  III,  para  32 

TMIN 

Minimum  wave  period  to  avoid  gen¬ 
eration  of  spurious  waves,  sec 

See  Part  III,  para  35 

KH 

Nondlmensional  water  depth 

HSR 

Height-to-stroke  ratio  for  two- 
dimensional  transfer  function 

See  Part  III,  para  21 

HO 

Wave  height  based  on  height-to- 
stroke  ratio,  in. 

HT 

Wave  height  based  on  three- 
dimensional  transfer  function 

See  Part  III,  para  23 

HMAX 

Maximum  prebreaking  wave  height 

See  Part  III,  para  24 

HLR 

Wave  steepness  ratio 

Measure  of  nonlinearity 
see  Part  III,  para  39 

PNL 

Goda's  dimensionless  nonlinear 
parameter 

Measure  of  nonlinearity 
see  Part  III,  para  39 

TIME1 

Time  prior  to  data  collection 

See  Part  III,  para  49 

TIME2 

Time  for  specified  wave  to  travel 

See  Part  III,  para  49 

to  back  wall  and  return 


I 


Table  8 

Summary  of  Input  Parameters  for  Process  IDCAL 


Header  2  Variables 


Line  1 
Period <  sec 

Line  2 
Scans /Period 

Line  3 

Records /Period 

Lin< 
Number  oi 

0.75 

38 

38 

6  i 

1.00 

50 

50 

41 

1.50 

75 

75 

1< 

2.00 

100 

100 

< 

3.00 

150 

150 

- 

Line  16 


Table  9 

Relative  Times  for  Data  Collection 


Number  of 

Test 

Generator 

Measurement 

Period 

Test 

Duration 

Start  Delay 

Delay 

T 

sec 

Periods 

sec 

sec 

sec 

Table  10 

Absolute  Times  for  Data  Collection 


Table  11 


Su—ary  of  Teat  Caae  Parameters.  Monochrogatic 
Perforaance  Testa  In  DSWG  Baeln 


Test  Duration 

Travel 

Tine 

Data 

Collection 

Teat 

ID 

Period 

sec 

Number 
of  vaves 

sec 

t 

g 

fpa 

Tlael 

see 

Ti«e2 

sec 

Nunber 

Begin 

sec 

End 

sec 

MPTST9 

0.75 

120 

90 

2.07 

14.5 

74.4 

64 

30 

78 

MPTST1 

1.00 

90 

90 

3.04 

9.9 

50.6 

48 

24 

72 

MPTST6 

1.50 

60 

90 

4.29 

7.0 

35.9 

16 

21 

45 

HPTST2 

2.00 

60 

120 

4.86 

6.2 

31.7 

6 

16 

28 

MPTST3 

3.00 

20 

60 

5.29 

5.7 

29.1 

3 

13 

22 

Table  12 

Wave  Period  Analysis 


Waves  of  0  deg 

MONOCHROMATIC  PERFORMANCE  TESTS  IN  DSWG  BASIN 
Measured  vs.  Theoretical  Wave  Periods 
Direction  *  0  Deg,  Depth  =  1  Ft 

February/March  86 

Overall  BRF, 

•/.: 

99.65 

S=1M  S=3" 

Period  Meas  Meas 

( Sec  >  ( Sec  >  ( Sec  > 

S=6“  Sz9" 
Meas  Meas 
(Sec)  (Sec) 

S=  1 1 " 

Meas 

(Sec) 

Meas 

Ave 

(Sec) 

BRF 

(*/.) 

0.75 

0.740 

0.740 

0.740 

98.67 

1.00 

1.000 

1.000 

1.000 

1.000 

100. 00 

1.50 

1.501 

1.500 

1.502 

1.502 

1.501 

99.92 

2.00 

1.999 

1.999 

1.999 

1.997 

1.998 

1.998 

99.92 

3.00 

2.974 

2.994 

3.002 

2.998 

2.997 

2.993 

99.77 

Waves  of  15  deg 

MONOCHROMATIC  PERFORMANCE  TESTS  IN  DSWG  BASIN 

Measured  vs 

.  Theoretical  Wave 

Periods 

Direction  = 

15  Deg, 

Depth  =  1 

Ft 

Overall  BRF 

99.61 

February/March  86 

S*l" 

S=3" 

S*6"  S=9” 

S=  1 1 " 

Meas 

Period 

Meas 

Meas 

Meas  Meas 

Meas 

Ave 

BRF 

(Sec) 

(Sec) 

(Sec ) 

(Sec)  (Sec) 

(Sec) 

(Sec) 

C/.) 

0.75 

0.740 

0.740 

0.740 

98.67 

1.00 

1.000 

1.000 

1.000 

100. 00 

1.50 

1.499 

1.500 

1.500 

1.500 

99.90 

2.00 

1.999 

1.999 

2.031  1.997 

2.007 

99.68 

3.00 

2.996 

2.982 

3.000  2.987 

2.994 

2.992 

99.  73 

(Continued) 
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Table  12.  (Continued) 


Waves  of 

30  deg 

MONOCHROMATIC  PERFORMANCE  TESTS  IN  DSWG  BASIN 

Measured  vs 

.  Theoretical  Wave 

Periods 

Direction  = 

30  Deg , 

Depth  =  1 

Ft 

Overall  BRF 

,  7.: 

99.63 

February/March  86 

S-l" 

S=3" 

S=6" 

S*9" 

S=ll" 

Meas 

Per i od 

Meas 

Meas 

Meas 

Meas 

Meas 

Ave 

BRF 

(Sec) 

(Sec) 

(Sec) 

(Sec) 

(Sec) 

(Sec) 

(Sec) 

(7.) 

0.75 

0.740 

0.740 

98.67 

1.00 

1.000 

1.000 

1.000 

100. 00 

1.50 

1.501 

1.500 

1.499 

1.500 

100. 00 

2.00 

2.002 

2.000 

2.028 

1.998 

2.007 

99.65 

3.00 

2.992 

3.003 

3.  003 

2.990 

2.988 

2.995 

99.84 

Waves  of  45  deg 

MONOCHROMATIC  PERFORMANCE  TESTS  IN  DSWG  BASIN 


Measured  vs.  Theoretical  Wave  Periods 


Direction  * 

45  Deg , 

Depth  = 

1  Ft 

February/March  86 

S*l" 

S=3" 

S=6" 

Per i od 

Meas 

Meas 

Meas 

(Sec) 

(Sec) 

(Sec) 

(Sec) 

0.75 

0.740 

1.00 

1.000 

0.998 

1.50 

1.500 

1.499 

1.495 

2.00 

1.999 

1.996 

1.992 

3.00 

2.993 

2.979 

2.977 

Overall  BRF ,  7.:  99.51 


S<=9" 

Meas 

(Sec) 

S=  1 1 " 
Meas 
(Sec) 

Meas 

Ave 

(Sec) 

BRF 

(7.) 

0.740 

98.67 

0.999 

99.90 

1.498 

99.87 

2.002  . 

1.997 

99.86 

2.978 

2.962 

2.978 

99.26 

(Continued) 
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Table  12.  (Concluded) 


Waves  of  60  deg 

MONOCHROMATIC  PERFORMANCE  TESTS  IN  DSWG  BASIN 


Measured  vs.  Theoretical  Wave 
Direction  =  60  Deg,  Depth  =  1 
February/March  86 

Per iods 
Ft 

Overall  BRF, 

•/.: 

99.33 

S=l" 

Period  Meas 

<  Sec )  ( Sec ) 

S=3" 

Meas 

(Sec) 

S=6" 

Meas 

(Sec) 

S=9” 

Meas 

(Sec) 

S=1 1 " 

Meas 

(Sec) 

Meas 

Ave 

(Sec) 

BRF 

(•/.) 

0.75 

0.740 

0.740 

98.67 

1.00 

1.000 

0.997 

0.999 

99.85 

1.50 

1.501 

1.497 

1.497 

1.498 

99.89 

2.00 

1.993 

1.989 

1.992 

1.985 

1.990 

99.49 

3.00 

2.999 

3.003 

2.990 

2.935 

2.883  2.962 

98.73 
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Table  13 

Measured  Average  Wave  Periods  for  All  Wave  Directions 


MONOCHROMATIC  PERFORMANCE  TESTS  IN  DSWG  BASIN 
Measured  Average  vs.  Theoretical  Wave  Periods 

Depth  =  1  Ft  Overall  BRF ,  '/.:  99.61 


February/March  B6 


Period 

(Sec) 

Dir=0 
Ave 
(Sec ) 

Dir=15 

Ave 

(Sec) 

Di  r=30 
Ave 
(Sec) 

0.75 

0.740 

0.740 

0.740 

1.00 

1.000 

1.000 

1.000 

1.50 

1.501 

1.500 

1.500 

2.00 

1.99S 

2.007 

2.007 

3.00 

2.993 

2.992 

2.995 

Di r =45 
Ave 
(Sec) 

Di  r=60 
Ave 
(Sec ) 

Overal 1 
Ave 
(Sec) 

BRF 

(•/.) 

0.740 

0.740 

0.740 

98.67 

0.999 

0.999 

1.000 

99.96 

1.498 

1 . 498 

1.499 

99.96 

1.997 

1.990 

2.000 

99.99 

2.978 

2.962 

2.984 

99.47 

* 


[ 


•  <, 


Table  14 

Harmonic  Analysis 


Waves  of  0  deg 


MONOCHROMATIC  PERFORMANCE  TESTS  IN  DSNS  8ASIN 
Har ionic  Analysis  of  Linear  Haves 
Percent  Variance  in  First  Harionic 
Direction  =  0  Deg  Depth  =  1  Ft 
February/April  86 


S=l* 

S=3* 

S=6* 

S=9* 

s=n* 

Period 

Soda 

Var 

Soda  Var 

Soda 

Var 

Soda 

Var 

Goda 

Var 

(Sec) 

NLP 

(X) 

NLP  (XI 

NLP 

(X) 

NLP 

(X) 

NLP 

(I) 

0.056 

99.6 

1.00 

0.035 

99.7 

0.105 

97.2 

1.50 

0.029 

99.8 

0.086 

98.5 

0.172 

88.0 

2.00 

0.031 

99.7 

0.093 

97.4 

0.187 

89.1 

0.280 

76.2 

0.342 

68.2 

3.00 

0.040 

99.7 

0.121 

96.7 

0.242 

83.9 

0.363 

68.5 

0.443 

57.9 

Waves  of  15  deg 

HQNOCHROHATIC  PERFORMANCE  TESTS  IN  DSN6  BASIN 
Harionic  Analysis  of  Linear  Haves 
Percent  Variance  in  First  Harionic 
Direction  =  IS  Deg  Depth  =  1  Ft 
February/April  86 


S=r  S=3*  S*6“  S=9* 


Period 

(Sec) 

6oda 

NLP 

Var 

(X) 

6oda 

NLP 

Var 

(1) 

6oda 

NLP 

Var 

\  %i 

Seda 

NLP 

Var 

(X) 

Soda 

NLP 

Var 

(I) 

0.75 

0.058 

99.6 

1.00 

0.036 

99.8 

0.109 

96.4 

1.50 

0.030 

99.9 

0.089 

98.7 

0.178 

87.8 

2.00 

0.032 

99.6 

0.097 

97.2 

0.193 

89.0 

0.290 

75.8 

3.00 

0.042 

99.4 

0.125 

96.5 

0.250 

84.4 

0.375 

68.4 

0.459 

56.4 

1 
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Table  14.  (Continued) 


Waves 

MONOCHROMATIC  PERFORMANCE  TESTS  IN  DSNS  BASIN 
Haraonic  Analysis  of  Linear  Naves 

Percent  Variance  in  First  Haraonic 

Direction  =  30  Deg  Depth  =  1  Ft 

February /April  86 

S=l*  S=3‘ 

of  30  deg 

S*6* 

S=9* 

S=ll* 

Period 

Soda 

Var 

6oda 

Var 

Soda 

Var 

Soda 

Var 

Soda 

Var 

(Sec) 

NLP 

(Z) 

NLP 

(Z) 

NLP 

(Z) 

NLP 

(Z) 

NLP 

(Z) 

0.75 

0.064 

99.6 

1.00 

0.041 

99.8 

0.122 

97.5 

1.50 

0.034 

99.9 

0.101 

98.4 

93.6 

2.00 

0.03b 

96.2 

0.109 

96.2 

0.217 

82.2 

0.326 

67.7 

3.00 

0.047 

99.1 

0.140 

96.3 

0.280 

82.6 

0.420 

64.2 

0.513 

52.9 

Waves 

of  45  deg 

MONOCHROMATIC  PERFORMANCE  TESTS  IN  0SN6  BASIN 

Haraonic  Analysis  of  Linear  Haves 

Percent  Variance  in  First  Haraonic 

Direction 

*  45  Deg 

Depth 

=  1  Ft 

February/April  66 

s=r 

S=3* 

S*6" 

S=9' 

S=ll* 

Period 

Soda 

Var 

Soda 

Var 

Soda 

Var 

Soda 

Var 

Soda 

Var 

(Sec) 

NLP 

(Z) 

NLP 

(Z) 

NLP 

!Z) 

NLP 

(Z) 

NLP 

(Z) 

0.75 

0.072 

99.4 

1.00 

0.053 

99.8 

0.128 

95.8 

1.50 

0.042 

99.6 

0.127 

96.4 

2.00 

0.045 

99.4 

0.134 

94.3 

0.268 

77.0 

3.00 

0.05 6 

99.1 

0.168 

92.8 

0.335 

70.4 

0.503 

36.8 

0.614 

36.4 

(Continued) 
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Table  14.  (Concluded) 


Waves  of  60  deg 


MONOCHROMATIC  PERFORMANCE  TESTS  IN  DSM6  BASIN 
Haraonic  Analysis  of  Linear  Haves 
Percent  Variance  in  First  Haraonic 
Direction  =  60  Deg  Depth  =  I  Ft 
February/April  86 


S=J'  S=6"  S=9‘ 


Period 

(Sec) 

Soda 

NLP 

Var 

(!) 

Soda 

NLP 

Var 

(X) 

6oda 

NLP 

Var 

(2) 

Soda 

NLP 

Var 

(Z) 

0.75 

0.072 

95.3 

1.00 

0.053 

99.3 

1.50 

0.056 

99.5 

0.169 

93.9 

2. 00 

0.070 

97.5 

0.211 

87.3 

0.362 

53.7 

3.00 

0.077 

98.8 

0.231 

91.3 

0.463 

64.6 

0.694 

40.9 
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Table  15 

Wave  Height  Analysis 


Waves  of  0  deg 

HOMOCHAOHhTIC  PERFORNANCE  TESTS  IN  DSNS  BASIN 
Neasured  vs.  Theoretical  Nave  Heights 

Direction  »  I  Deg,  Depth  *  I  Ft  Overall  BRF,  I:  87.45 

February  84 


Period 

(Sec) 

Nax 

Height 

(In) 

Stroke  *  1* 
Neas  Theory 

(In)  (In) 

Stroke  •  3’ 
Neas  Theory 
(In)  (In) 

Stroke  *  4* 
Reas  Theory 

(In)  1  In) 

Stroke  *  9* 
Neas  Theory 
(In)  (In) 

Stroke  ‘11* 
Neas  Theory 

(In)  (In) 

B.75 

3.31 

1.55 

1.77 

3.88 

3.38 

l.BB 

4.79 

1.17 

1.31 

3.27 

3.94 

i.SI 

4.23 

B.4B 

8,81 

1.77 

2.42 

3.95 

4.83 

2.BB 

4.79 

1.41 

B.58 

1.33 

1.75 

2.88 

3.49 

5.14 

5.24 

4.38 

4.48 

3.BB 

7.2B 

B.23 

B.3B 

B.B1 

1.13 

1.91 

2.27 

3.22 

3.48 

4.24 

4.14 

BRF,  Z: 

79.59 

81.82 

B1.78 

95.83 

99.24 

Notes: 

1.  Slight  Breaking  Occurred 

2.  BD  *  Bad  Data 


Waves  of  15  deg 


hGNOChRGNAT iC  PERFORKftHCE  TESTS  IN  DSH6  BASIN 
Neasured  vs.  Theoretical  Nave  Heights 

Direction  =  15  Deg,  Depth  =  1  Ft  Overall  BRF,  I:  8a. 48 

Fehruary/Harch  84 


Nax  Stroke  »  1’  Stroke  =  3'  Stroke  =  4’  Stroke  =  V  Stroke  =  U1 


Period 

(Sec) 

Height 

(In) 

Neas 

(In) 

Theory 

(In) 

Neas 

(In) 

Theory 

(In) 

Neas 

(In) 

Theory 

(In) 

Neas 

(In) 

Theory 

(In) 

Neas 

(In) 

Theory 

(In) 

1.75 

3.38 

1.4B 

1. 84 

l.BB 

4.79 

1.84 

1.34 

3.34 

4.87 

1.58 

4,23 

8.54 

8.83 

1.83 

2.58 

4.24 

5.88 

2.88 

4.79 

8.38 

8.48 

1.38 

1.81 

2.87 

3.41 

5.88 

5.42 

3.18 

7.28 

8.25 

8.39 

8.85 

1.17 

1.94 

2.35 

3.37 

3.52 

4.44 

4.38 

BRF,  I: 

73.98 

74.45 

82.94 

93.42 

94.28 

Notes: 

1.  Slight  Breaking  Occurred 

2.  BD  >  Bad  Data 


(Continued) 
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Table  IS.  (Continued) 


Waves  of  30  deg 


HGhOCHAOHATIC  PERFGRHANCE  TESTS  IN  D5H6  BASIN 
Heasured  vs.  Theoretical  Nave  Heights 
Direction  -  31  Be),  Depth  *  t  Ft 
February/Harch  84 


Ha*  Stroke  *  1*  Stroke  •  3* 


Period 

(Sec) 

Height 

(In) 

Neas 

(In) 

Theory 

(In) 

Heas 

(In) 

Theory 

(In) 

1.75 

3.31 

1.34 

2.B1 

l.BI 

4.79 

1.11 

1.52 

3.3B 

4.55 

l.SI 

4.23 

B.59 

B.94 

1.9B 

2.B2 

2.BI 

4.79 

B.4B 

B.48 

1.4B 

2,13 

3.H 

7.2B 

B.2B 

B.44 

B.B2 

1.31 

IF,  I: 

45.47 

7B.B3 

Overall  BRF,  Z: 

B1.4B 

Stroke  •  4* 

Stroke  *  9* 

Stroke  *11’ 

Heas  Theory 

Heas  Theory 

Heas  Theory 

(In)  (In) 

(In)  1  In) 

(In)  1 In) 

4.49 

5.44 

3.42 

4.B7 

5.42 

4.11 

2.BB 

2.42 

3.45 

3.94 

4.BB 

4.81 

81.37 

92.33 

99.79 

Note*: 

1.  Slight  Breaking  Occurred 

2.  BO  *  Bad  Data 


Waves  of  45  deg 


NONOCHRONATIC  PERFORHANCE  TESTS  IN  DSNS  BASIN 
Heasured  vs.  Theoretical  Nave  Heights 

Direction  *  45  Deg,  Depth  =  1  Ft  Overall  BRF,  Z:  49.21 

February/Harch  84 


Ha*  Stroke  *  1*  Stroke  =  3*  Stroke  *  6'  Stroke  *  9'  Stroke  *  11* 


Period 

(Sec) 

Height 

lint 

Heas 

(In) 

Theory 

(In) 

Heas 

(In) 

Theory 

(In) 

Heas 

(In) 

Theory 

(In) 

Heas 

(In) 

Theory 

(In) 

Neas 

(In) 

Theory 

(In) 

1.75 

3.31 

1.32 

2.27 

l.BB 

4.79 

1.15 

2.N 

4.14 

4.79 

1.51 

4.23 

1.71 

1.19 

2.51 

3.57 

5.13 

4.23 

2. IB 

4.79 

1.41 

1.84 

1.55 

2.51 

4.14 

5.12 

4.54 

4.79 

3.11 

7.21 

1.24 

1.52 

1.93 

1.57 

2.53 

3.14 

BD 

4.71 

BD 

5.74 

BRF,  Z: 

54.11 

72.51 

81.45 

44.84 

Notes: 

1.  Slight  Breaking  Occurred 

2.  BD  •  Bad  Data 
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(Sheet  2  of  3) 


T 


Table  15.  (Concluded) 


Waves  of  60  deg 

HONOCHflOHATlC  PERFORHANCE  TESTS  IN  DSNS  BASIN 
Heasured  »*•  Theoretical  Nave  Heights 

Direction  *  61  Deg,  Depth  *  1  Ft  Over  ill  Biff,  I:  72.71 

Fehruary/Harch  86 


Hu  Stroke  •  1*  Stroke  •  3*  Stroke  •  6'  Stroke  •  9*  Stroke  •  11* 


Period 

(Sec) 

Height 

(In) 

Hens 

(In) 

Theory 

(In) 

Heis 

(In) 

Theory 

(In) 

Kees 

(In) 

Theory 

(In) 

fleas 

(In) 

Theory 

(In) 

Heas 

(In) 

Theory 

lie) 

I.7S 

3.31 

1.92 

2.27 

1.11 

4.7? 

1.16 

2.M 

4.51 

4.79 

l.SI 

6.23 

1.92 

1.58 

3.17 

4.73 

2.11 

6. 79 

I.S4 

1.32 

2.22 

3.95 

4.58 

6.79 

3.11 

7.21 

1.28 

8.72 

1.37 

2.17 

3.16 

4.34 

5.37 

6.51 

5.71 

7.21 

BRF,  X: 

61.19 

71.99 

68.82 

B2.49 

79.17 

Note: 

1.  Slight  Breaking  Occurred 

2.  ID  *  lid  Data 
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Table  16 

Breaking  Wave  Tests 


MCNCCHROHATIC  PERFORMANCE  TESTS  IN  DSNS  BASIN 
Breaking  Have  Tests 
Depth  =  1  Ft 
February  86 


Theory  Beloe  Breaking  Near  Breaking  Above  Breaking 


Period 

(Sec) 

Stroke 

(In) 

Height 

(In) 

Stroke 

(In) 

Height 

(In) 

Stroke 

(In) 

Height 

(In) 

Stroke 

(In) 

Height 

(In) 

8.75 

1.86 

3.38 

2.85 

2.44 

2.41 

2.95 

2.78 

3.81 

1.88 

3.45 

4.79 

4.13 

4.45 

4.39 

4.75 

4.68 

4.95 

1.58 

7.4? 

4.23 

7.49 

5.35 

8.46 

6.14 

8.78 

6.18 

2.88 

11.78 

4.79 

18.21 

5.74 

11.28 

6.38 

11.71 

6.37 

3.81 

18.91 

7.28 

11.91 

4.67 

Table  17 

Constancy  of  Measured  Wave  Heights 


Measured  Nave  Heifkt,  In 


Theory 

l  •  11  Ft 

I  *  21  Ft 

I  «  38  Ft 

AH 

(X) 

(Sk) 

(Dog) 

ntiyii 

(la) 

33  43  S3 

33  43  53 

33 

45  33 

ii 

l.N 

1 

3.94 

2.932 

2.953 

3.873 

Stall  Nave  Safes 

3.394  3.338  3.442 

3.243 

3.213 

3.829 

3.314 

7.8 

IS 

4.17 

2.732 

2.491 

3.334 

3.383 

3.334 

3.491 

3.822 

3.299 

3.724 

3.348 

8.8 

31 

4.33 

2.924 

3.179 

3.894 

3.328 

3.372 

2.984 

3.983 

3.881 

3.343 

3.249 

7.9 

45 

Bfi 

4.224 

4.348 

3.943 

4.185 

3.748 

3.348 

3.418 

2.888 

2.348 

3.424 

14.4 

1.31 

• 

2.42 

2.828 

1.774 

1.778 

1.433 

1.371 

1.793 

1.971 

1.941 

1.939 

1.824 

7.1 

IS 

2.31 

1.824 

2.881 

1.934 

1.712 

2.893 

1.827 

2.182 

1.599 

1.755 

1.881 

8.1 

31 

2.B2 

2.157 

1.883 

2.898 

2.883 

2.834 

2.184 

1.719 

1.998 

2.218 

2.819 

3.8 

43 

3.37 

2.244 

2.474 

2.917 

3.244 

2.849 

2.844 

2.259 

1.948 

1.528 

2.415 

18.4 

Jordan  Nave  6a|ti 


l.N 

8 

3.94 

3.138 

3.9N 

2.343 

3.338 

3.439 

3.281 

3.344 

2.735 

3.934 

3.236 

9.7 

IS 

4.87 

3.813 

3.747 

3.383 

3.183 

3.357 

3.288 

2.848 

3.234 

3.934 

3.331 

7.3 

38 

4.53 

2.791 

3.477 

2.744 

3.8N 

3.778 

2.943 

4.843 

2.778 

3.339 

3.382 

13.1 

43 

BR 

4.313 

4.477 

4.244 

3.724 

3.886 

3.137 

3.434 

2.416 

2.216 

3.579 

17.2 

1.38 

a 

2.42 

1.828 

1.481 

1.728 

1.775 

1.387 

1.784 

1.723 

1.715 

1.739 

1.712 

3.1 

13 

2.58 

1.427 

1.743 

1.773 

1.787 

2.172 

1.822 

1.883 

1.514 

1.628 

1.743 

6.7 

38 

2.82 

1.983 

1.812 

1.943 

2.823 

1.944 

2.123 

1.493 

1.934 

2.844 

1.926 

4.3 

43 

3.57 

2.831 

2.513 

2.745 

3.844 

2.712 

1.981 

1.932 

1.632 

1.238 

2.287 

22.2 

Notes: 

1.  M  *  Nave  height  should  exceed  saaiouo  possikli  pre-breakiaj  uavi  beifht. 


Table  18 

Wave  Direction  Analysis 


' 

( 

Waves  of  0  deg 

MDMOCHROHATIC  PERFORMANCE  TESTS  IN  DSNS  BASIN 
Measured  vs.  Theoretical  Have  Directions 


Direction  =  0  Deg,  Depth  1  1  Ft 

Overall  BRF,  l: 

99.9 

February/March  86 

Theory 

s=r 

S-3’ 

Ss6* 

3=9" 

5=11" 

Meas 

Period  Direction 

Meas 

Meas 

Meas 

Meas 

Meas 

Ave 

BRF 

(Sec)  (Deg) 

(Deg) 

(Deg) 

(Deg) 

(Deg) 

(Deg) 

(Deg) 

(l) 

8.75 

8.8 

8.8 

8.8 

8.8 

mm 

1.88 

8.8 

8.8 

8.7 

8.8 

8.2 

HE 

1.58 

8.8 

8.8 

8.8 

8.8 

8.8 

El 

2.88 

8.3 

8.6 

8.8 

8.3 

8.6 

8.4 

99.9 

5.88 

8.8 

8.6 

i.3 

8.8 

8.8 

1.3 

e.B 

99.8 

Waves  of  15  deg 

HONOCHRGMATIC  PERFORMANCE  TESTS  IN  DSHG  BASIN 
Measured  vs.  Theoretical  Have  Directions 


Direction  =  15  Deg,  Depth  =  1  Ft 

Overall  BRF,  Is 

97.9 

February/Narch  B6 

Theory 
Period  Direction 
(Sec)  (Deg) 

S=i" 

Meas 

(Deg) 

S=3" 

Meas 

(Deg) 

S*6" 

Meas 

(Deg) 

S=9" 

Meas 

(Deg) 

S=ll" 

Meas 

(Deg) 

Meas 

Ave 

(Deg) 

BRF 

(1) 

8.75 

15.6 

15.4 

15.4 

15.4 

98.7 

1.18 

14.5 

>4.6 

14.6 

14.6 

99.3 

1.58 

14.9 

15.9 

16.8 

15.9 

15.9 

93.1 

2.88 

14.9 

16.5 

13.2 

15.2 

15.2 

15.8 

99.2 

3.88 

14.9 

14.3 

15.3 

15.8 

15.5 

15.8 

99.2 

(Continued) 
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Table  18.  (Continued) 


Waves  of  30  deg 

MONOCHROMATIC  PERFORMANCE  TESTS  IN  0SM6  BASIN 
Measured  vs.  Theoretical  Nave  Directions 


Direction  =  38  Deg,  Depth  =  1  Ft 

Overall  BRF,  I; 

97.2 

February/March  B6 

Theory 
Period  Direction 
(Seel  (Deg) 

s=r 

Heas 

(Deg) 

S=3* 

Heas 

(Deg) 

S=6“ 

Meas 

(Deg) 

S=9’ 

Heas 

(Deg) 

S=ll" 

Meas 

(Deg) 

Meas 

Ave 

(Deg) 

BRF 

(X) 

0.75 

28.8 

29.4 

29.4 

94.3 

1.88 

38.2 

38.2 

31. 8 

31.8 

97.4 

1.58 

31.8 

31.5 

38.7 

29.9 

38.7 

99.8 

2.88 

38.9 

32.3 

38.5 

33.8 

32.3 

32.8 

94.4 

3.88 

38.3 

34.1 

28.9 

29.9 

38.8 

38.4 

38.7 

9B.8 

Waves  of  45  deg 


MCNOCHROMATIC  PERFORMANCE  TESTS  IN  DSNS  BASIN 
Measured  vs.  Theoretical  Nave  Directions 

Direction  =  45  Deg,  Depth  =  I  Ft  Overall  BRF,  X:  9B.8 

February/March  86 


Theory 

s=r 

5=3’ 

5=6’ 

S=9* 

s=ir 

Meas 

Period  Direction 

Meas 

Meas 

Meas 

Meas 

Meas 

Ave 

BRF 

(Sec)  (Deg) 

(Deg) 

(Deg) 

(Deg) 

IDeg) 

(Deg) 

(Deg) 

(X) 

0.75 

38.8 

38.8 

38.8 

188.8 

1.88 

48.9 

51.1 

4B.2 

49.7 

98.5 

1.58 

47.4 

47.0 

47.6 

4B.3 

47.6 

99.5 

2.88 

45.9 

44.2 

44.2 

46.6 

45.0 

98.8 

3.81 

43.9 

43.7 

44.1 

45.1 

45.9 

44.6  44.7 

98.2 
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Table  18.  (Concluded) 


Waves  of  60  deg 


MONOCHROMATIC  PERFORMANCE  TESTS  IN  DSN6  BASIN 
Measured  vs.  Theoretical  Have  Directions 


Direction  *  61  Deg,  Depth  •  1  Ft 

Overall  BRF,  I: 

98.4 

February/March  86 

Theory 
Period  Direction 
(Sec)  (Deg) 

S=l* 

Meas 

(Deg) 

S=3* 

Meas 

(Deg) 

S*6* 

Meas 

(Deg) 

Ss9’ 

Meas 

(Deg) 

s=ir 

Meas 

(Deg) 

Meas 

A  ve 
(Deg) 

BRF 

(Z) 

8.75 

38.8 

38.8 

38.8 

188.8 

1.88 

48.9 

49.1 

49.1 

49.1 

99.6 

1.58 

59.2 

59.6 

59.7 

59.7 

99.2 

2.88 

65.8 

62.8 

62.8 

62.4 

97.8 

3.00 

58.5 

54.5 

57.3 

58.8 

52.8 

55.7 

95.1 
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Table  19 

Measured  Average  Versus  Theoretical  Wave  Direction  for  All  Directions 


mOCWtOMTIC  PERFQRHANCE  TESTS  IN  DSNS  BASIN 
hiiirS  Average  vs.  Theoretical  Navi  Directions 

Beyth  *  1  Ft  Ov trail  BRF,  2:  96.81 

Felraary/Harch  86 


Perio* 

(See) 

Hu  Birtction 
Anflt  1  Reg'd 
(Deg)  PuMltt 

Direction  *  8 
Heat  Theory 
(Dig)  (Beg) 

Direction  *  IS 
Heu  Theory 
(Deg)  (Deg) 

Directioi  ■  38 
Heu  Theory 
(Deg)  (Deg) 

Direction  *  43 
Hens  Theory 
(Deg)  (Deg) 

Direction  *  68 
Hens  Theory 
(Deg!  (Deg) 

6.73 

38.8 

3 

8.6 

8.8 

13.4 

15.6 

29.6 

28.6 

38.  B 

38.8 

38.8 

38.8 

l.N 

48.9 

4 

8.2 

8.8 

14.6 

14.3 

31.8 

38.2 

49.7 

48.9 

49.1 

48.9 

1.31 

39.2 

6 

8.8 

6.8 

13.9 

14.9 

38.7 

31.8 

47.6 

47.4 

59.7 

59.2 

2.M 

63.8 

8 

8.4 

8.6 

15.8 

14.9 

32.6 

38.9 

45.6 

45.9 

62.4 

63.8 

3.M 

67.5 

12 

6.8 

8.8 

IS.  8 

14.9 

38.7 

36.3 

44.7 

43.9 

55.7 

58.5 

BRF,  X: 

99.6 

97.6 

99.6 

98.7 

Table  20 

Wave  Direction  Verification  Tests 


MONOCHROMATIC  PERFORMANCE  TESTS  IN  DSWG  BASIN 
Direction  Tests:  Measured  vs.  Theoretical  Wave  Direction 
Desired  Direction,  60  Deg 
February/March  36 


T  heur  y 

5”  3" 

Per i o J 

Di recti  on 

Me  as 

DRF 

(Sfc  > 

( Dog  > 

<  Deg  > 

(■/.:• 

1.50 

50 . 3 

52.9 

94. 8 

1.50 

59.2 

59.6 

99.3 

1 . 50 

69.6 

70.2 

99.  1 

Overall  BRF,  % 


97. C 


Table  21 

Nonlinear  Wave  Period  Analysis 


Wave 8  of  0  deg 

MONOCHROMATIC  PERFORMANCE  TESTS  IN  DSWG  BASIN 
Nonlinear  Measured  vs.  Theoretical  Have  Periods 


Direction  =  0  Deg, 
February /March  86 

Depth  = 

1  Ft 

Overall  BRF, 

X: 

99.67 

Period 

(Sec) 

S=6" 

Meas 

(Sec) 

S=9" 

Meas 

(Sec) 

S=ll“ 

Meas 

(See) 

Meas 

Ave 

(Sec) 

BRF 

(*) 

1.50 

2.00 

3.00 

1.498 

1.993 

3.017 

1.495 

1.994 

3.016 

1.500 

1.997 

1.498 

1.995 

3.017 

99.84 

99.73 

99.43 

Waves  of  30  deg 

MONOCHROMATIC.  FERFORMANCE  TESTS  IN  BSWG  BASIN 

Nonlinear  Measured  vs.  Theoretical  Wave  Periods 


Direction  =  30  Deg, 
February/March  86 

Depth  = 

1  Ft 

Overall  BRF, 

7.: 

99.71 

Per i od 
(Sec ) 

S-6” 

Meas 

(Sac) 

S-9" 

Meas 

(Sac) 

S*  1 1  ” 

Meas 

(Sac) 

Meas 

Ave 

(Sac) 

BRF 

(7.) 

1.50 

2.00 

3.00 

1.497 

1.997 

2.999 

1.497 

1.989 

2.999 

1 . 494 
1.979 

1.496 

1 .938 
2.999 

99.73 
99.  42 
99.97 

Table  22 

Comparison  of  Linear  and  Nonlinear  Harmonic  Analysis  Results 


MONOCHROMATIC  PERFORMANCE  TESTS  IN  DSWG  BASIN 

Comparison  o-f  Harmonic  Analysis 

Linear  vs  Nonlinear  Waves 

Depth  -  1  Ft 

February/April  96 


D=0 

Goda 

NLP 

Mono 

Nonl i near 

0.  172 

88.  0 

87.  1 

0.  187 

88.  1 

86.7 

0.  242 

03.9 

75.  4 

0.  280 

76.  2 

78.5 

0.  342 

68.  2 

74.4 

0.  363 

68.  5 

69.9 

D--30 

Goda 

NLP  Mono  Nonlinear 


0.201  93.6  91.7 

0.217  82. 2  88. 

0.280  82.6  79. 

0.326  67.7  79. 

0.420  64.2  71.1 


Nl  <3  Nl 


Table  23 

Nonlinear  Wave  Height  Analysis 


MONOCHROMA’: 
Nonlinear  1 
Direction  : 
February /Mt 

Period 

(Sec) 

Waves  of  0  deg 

[■IC  PERFORMANCE  TESTS  IN  DSWG  BASIN 

Measured  vs.  Theoretical  Wave  Heights 
:  0  Deg,  Depth  =  1  Ft 

irch  86  Overall  BRF,  X: 

S=6  S=9  S 

52.82 

=  11 

Meas 

(In) 

BRF 

(*) 

Meas 

(In) 

BRF 

(*) 

Meas 

(In) 

BRF 

(*) 

1.50 

3.63 

60.50 

5.13 

57.00 

6.15 

55.91 

2.00 

3.23 

53.83 

4.79 

53.22 

5.66 

51.45 

3.00 

2.63 

43.89 

4.13 

45.86 

BRF,  X: 

52.74 

52.03 

53.68 

Notes : 

1 .  Slight  Breaking 

2.  All  Breaking 

Waves 

of  30  deg 

MONOCHROMATIC  PERFORMANCE  TESTS  IN  DSWG  BASIN 

Nonlinear  Measured  vs 

.  Theoretical  Wave  Heights 

Direction  = 

=  30  Deg,  Depth  =  1 

Ft 

February/March  86 

Overall  BRF,  X: 

54.40 

S=6 

S 

=9 

S 

=  11 

Period 

Meas 

BRF 

Meas 

BRF 

Meas 

BRF 

(Sec) 

(In) 

(7.) 

(In) 

(X) 

(In) 

(X) 

1.50 

3.86 

64.38 

5.83 

64.76 

5.41 

49.  18 

2.00 

3.34 

55.68 

4.94 

54.93 

6.07 

55.  15 

3.00 

2.  76 

46.00 

4.  26 

47.36 

BRF,  V.s 

55.36 

55.68 

52.  17 

Notest 

1.  Slight  Breaking 

2.  All  Breaking 

Q  L>. 


Table  24 

Nonlinear  Wave  Direction  Analysis 


Waves  of  0  deg 

MONOCHROMATIC  PERFORMANCE  TESTS  IN  DSWG  BASIN 
Nonlinear  Measured  vs.  Theoretical  Wave  Directions 
i recti  on  -  0  Deg ,  Depth  ~  1  Ft 

ebruary /March  36  Overall  BRF ,  "/.  99.9 


Theory 

S~6" 

S--9" 

S—  1  i  " 

Meas 

Per i od 

Di recti  on 

Meas 

Meas 

Meas 

Ave 

BRF 

( Sec  > 

(Deg) 

(Deg ) 

(Deg  > 

( Deg ) 

( Deg ) 

('/.) 

1 . 50 

0.  0 

0.0 

0.  0 

100.0 

2.00 

0.0 

0.0 

0.  0 

0.0 

100.0 

3.00 

0.  0 

0.0 

0.8 

0.  9 

0.  6 

99.8 

Waves 

of  30  deg 

MONOCHROMA 

iTIC  PERFORMANCE  TEST 

S  IN  DSWG 

BASIN 

Nani  inear 

Measured 

vs.  Theoretical  Wave 

Di  recti  on 

3 

Di  r  ecti  on 

-  30  Deg , 

Depth  =  1 

Ft 

February /March  06 

Ovei 

rail  BRF, 

7.: 

Theory 

S-6  " 

b  —  V 

Meas 

Period  Direction 

Meas 

Meas 

Ave 

BRF 

( Sec ) 

(Deg ) 

(Deg) 

(Deg ) 

( Deg ) 

;■/.) 

1 . 50 

31 . 0 

34.0 

34 . 0 

90.  3 

2.00 

30.9 

35.  6 

43.2 

39.  4 

72.5 

3.00 

30.3 

33.2 

35.  9 

34.6 

86.0 

Photo  1.  Directional  spectral  wave  generator  basin 


j 


Photo  3.  Beach  wave  absorbers  opposite  the  DSWG 


Photo  7 


9.  Power  and  signal  conditioning 
controllers  and  amplifiers 


Photo  10.  Close-op  of  Getty  amplifiers 
and  transformers 


Photo  11.  Control  and  feedback  transducers 


Front  viav 


b.  Roar  viaw 

Photo  12.  Wava  ganarator  control  conaola 


Photo  13.  Syatea  control  panel 


Photo  14.  Servo 
aodula  panel 


Overhead  viev 


b.  Closeup  viev 
Photo  IS.  Have  |i|e  support  fri 
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....  J  iSSS«S**iRB2:3R 

Jill  I!——— 

SAKS  *  I  R R R » * X  9 A A 8 K X R * X 

*  1  ••••••#•••••••• 


RBXXAR=X2fiBX8RA 


i  1  ? 
,Ji*  ,  . 
111!  5 


44**44*4**44*'e4 


isistsst8sstt:« 


39X88S32S8SX2RB 


8X8  = 
Is  f  •  W 

i  X  R  8  S 
!  i***' 


KSS2XS38==S3S2A 


S8sa 

•f 


RSS8=83X3£8S888 


Ids: 
T  T  • 


ssssxesASSASssa 


8R28 

¥r-* 


38*8X8=82=88232 


8X2  = 
ft** 


R8=*»=2888888== 

*4444444*4*4**+ 


8  ARB 
TV* 


8888R888S833882 


8R23 

44** 


=8S2RS8S2R98s89 


*  A  a  3 

4  f  •  •’ 


8R23XXS83XS3R3S 
f TTTTTT *tTTTf Tf 


2X88 

ft** 


=2S*82====2S**R 


8X888SXRS8R8S8X 


8  R  X  2 

f  Is  •  • 


=38828888=388X8 
•  •  •  W  W  ^  ^  ^  w  w  •  • 


«  =  *2 

*  T  - 


£R=X3RS8S=S82SR 


88S8XX3=83S88SX  22X2 

|  *4  •  * 


E=S82SRR*3R8z3= I 3R8B 

8  W  X  8  X  •  *  -  f  W  y  f  f  W  X  (  *  i  •  • 


235*35352*^22*5 


X  R  8  X 

*4** 


"  83SX2SRXX8RX2XX 


3XXR 


=BSkEX8X=S83R38 


298  2 

*4** 


i  i 
liii 


A3 


Command  File  for  Program 
MONOSUMMARY 


*  SET  VERIFY 

*  ON  ERROR  THEN  GOTO  EXIT 

*  R  MONOSUMMARY1 
2 

5 

0  .  ,15a  , 30  a  ,45  a  , 60  a 
5 

la ,3. ,6. ,9. ,11a 
5 

a  75 ,1a  ,  1  a  5 , 2  a  ,3a 
1. 

30 

0 
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APPENDIX  C:  CONTROL  SIGNAL  GENERATION 


4 


Co— pd  Fll«  for  COMPONM 


•ON  CONTROL  Y  THEN  GOTO  EOF IT 
•DEFINE/USER  MODE  SYS* INPUT  SYS*COmAND: 

*R  [ OLOMCCLEAUE ] COMPONM 

•EOF IT i  DEF I NE/USER_MODE  SYS* INPUT  SYS* COMMAND: 
*  R  TWOEOF 

•ON  ERROR  THEN  GOTO  D1S2 
•DISH  MSAOt 
♦DIS2:  DISH  MSBO « 


Co—  nd  Flic  for  COMPONC4 


♦ON  CONTROL_Y  THEN  GOTO  EOF IT 

*DEF I NE/USER_MODE  SYS*INPUT  SYS*C0Mh^D: 

♦R  t  OLDMCCLEAUE  J  C0MP0NC4 

♦EOFITj  DEFINE/USER_MODE  SYS*INPUT  SYS*COMWNDj 
*  R  TWOEOF 

♦ON  ERROR  THEN  GOTO  DIS2 
♦DISH  MSAO: 

*DIS2i  DISM  MSBO: 


C3 


•  «  ' 


Sinusoidal  Control  Signal  Ganaratlon 


Log  on 
Osar  ID: 
Password : 

Changa  dlractory 


Actlvata  cc 
procedure 

Tapa  drive 


Hampton 


Gat  currant  password  from  Hsapton 

Sat  dafsult  to  [OLDMCCLEAVEJ 
(Nota:  Procadura  usad  during  thasc 
tests,  since  changed) 

Make  sura  taps  not  aountad  as  prograa 
mounts 

O-MSAO:,  1-MSBO: 


Control  signal  H17360 

filename 


Must  ba  6  characters  long 


Old  or  new  tapa 


Old  parameter  file  Y 


Changa  #  components  S 


Changa  periods 


OLD-Fllas  already  on  tape,  searches  for 
double  end-of-file  (EOF)  and  removes  In 
preparation  for  writing  new  file 

NEW-No  files  currently  on  tape,  writes 
new  file  @  beginning 

Y-Reads  from  FOR013.DAT  information  iron 
previous  use 

N-Use  If  no  previous  file  exists  ot  pre¬ 
viously  used  for  C0MP0NC4 

N-Uses  currant  vaJ</e  for  #  components 

Y>Querles  user  for  #  of  components  to  be 
combined 

N-Leavss  unchanged 

Y -Qua rise  user  for  periods  for  each 
component  specified,  can  be  Identical 
among  components 


(Continued) 


C4 


Table 


11  Change  peak-to-  Y 

peak  displacement 

Displacement  6.0 

12  Change  paddle  N 

offset  angle 

13  Parameters  menu  L 

14  Parameters  menu  Q 

15  Write  EOF's  Y 


16  Log  off 


LO 


(Concluded) 


_ Comments _ 

N-Leaves  unchanged 

Y-Querles  user  for  phases  for  each  com¬ 
ponent  specified,  usually  aero 

Y -Queries  user  for  stroke.  In. 

N-Leaves  unchanged 

N-Leaves  unchanged 

Y-Querles  user  for  offset  angle  neces¬ 
sary  to  generate  vave  angle 

L-Lists  specifications  for  each 
component 

Q-Beglns  writing  control  signal  to  tape, 
usually  run  2  min.  Enter  CTRL  C  to 
terminate 

Y-Writes  2  EOF’s  @  end  of  tape 
N-Use  only  if  filename  not  specified 


Tabic  C2 


Procedure  for  Program  C0MP0NC4 


Step  Description  Input  _  Comenta 


1  Log  on 

User  ID:  Hasp ton 

Password:  Get  current  password  from  Hampton 

2  Change  directory  SDOM  Set  default  to  [ OLDMCCLEAVE ] 

(Note:  Procedure  used  during  these 
tests,  since  changed) 

3  Activate  command  @NRO  Make  sure  tape  not  mounted  as  program 

procedure  mounts 

4  Tape  drive  0  O-MSAO:,  1-MSBO : 


5  Control  signal  MN7360  Must  be  6  characters  long 

filename 

6  Old  or  new  tape  Old  OLD-Files  already  on  tape,  searches  for 

double  end-of-file  and  removes  in  prep¬ 
aration  for  writing  new  file 

NEW«No  files  currently  on  tape,  writes 
new  file  @  beginning 

7  Old  parameter  file  Y  Y-Reads  from  FOROI3.DAT  information  from 

previous  use 

N-Use  if  no  previous  file  exists  or  pre¬ 
viously  used  for  COMPONM 

8  Change  #  components  N  N-Uses  current  value  for  #  components 

Y-Queries  user  for  #  of  components  to  be 
combined 

9  Depth  1.0  Enter  value  in  feet 

10  Change  wavelengths  Y  Y-Queries  user  for  wavelength  for  each 

component,  input  described  in  step  10a 
below 

N-Leaves  unchanged 

10a  Enter  wavelength  12.45  Corresponds  to  2.0  sec 

(Continued) 

C6 


Table  C2  (Concluded) 


Step 

Description 

Input 

Comments 

11 

Change  phases 

N 

N-Leavea  unchanged 

Y-Queries  user  for  phases  for  each  com¬ 
ponent  specified,  usually  zero 

12 

Change  peak-to- 
peak  displacement 

V 

Y-Querles  user  for  stroke.  In. 

N-Leaves  unchanged 

Displacement 

6.0 

13 

Change  paddle 

N 

N-Leaves  unchanged 

offset  angle 

Y*»Queries  user  for  offset  angle  neces¬ 
sary  to  generate  wave  angle 

14 

Parameters  menu 

L 

L«Lists  specifications  for  each 
component 

15 

Parameters  menu 

Q 

Q*Begins  writing  control  signal  to  tape, 
usually  run  2  min.  Enter  CTRL  C  to 
terminate 

16 

WTlte  EOF' s 

Y 

Y-Vrites  2  EOF' s  @  end  of  tape 

N-Use  only  if  filename  not  specified 

17 

Log  off 

LO 

C7 


Procedure  for 


Step  Description  Input 

1  Log  on 

User  ID:  Hampton 

Password: 

2  Type  input  file  T 

3  Input  filename  Ml/360 

4  Output  filename  Ml/360 

5  Main  menu  1 

6  Header  1  menu 

3:  Run  #  1 

4:  Model  name  M17360 

11 :  Wave  height  1.0 

7  Header  2  menu 

1:  Sampling  period  2 

2:  Scans/period  100 

3:  Records /period  100 

4:  it  Periods  6 

13:  DSWG  stroke  6 

16:  DAC  updates/per  40 

18:  Wave  direction  30 

8  Main  menu  7 

9  Main  menu  14 


10 


Main  menu 


15 


Le  D1 

)cess  IDCAL,  Setup 


Comments 


Get  current  password  from  Hampton 

T=Test,  M=Master 

Must  be  6  characters  long 

Saves  changes  to  same  file  input,  new 
file  if  given  different  name 

l=Header  1,  choose  from  1  to  15, 
repeats  when  exit  a  Header  submenu 

Select  line  corresponding  to  input 
desired,  repeats  after  each  entry, 
enter  -1  to  return  to  main  menu 


Should  match  step  4  above,  inches 

Select  line  corresponding  to  input 
desired,  repeats  after  each  entry, 
enter  -1  to  return  to  main  menu 

Sec 

Period  *  50 
Period  *  50 

Test  duration  =  #1  *  #4  -  12 
Same  as  Header  1  for  wave  height 
Period  *  20  (Note:  Actually  used  200, 
but  only  good  for  header,  so  did  not 
matter) 

Enter  0,  15,  30,  45,  or  60  deg 

Lists  rod  coefficients,  switches  should 
prevent  gage  rod  damage,  enter  -1  to 
return  to  main  menu 

Prints  all  data  entries  including  gage 
and  rod  coefficients,  enter  -1  to 
return  to  main  menu 


Exits  process  and  saves  to  file  speci¬ 
fied  as  output 


Table  D2 

Procedure  for  Proceaa  IDCAL,  Calibrating  Gages 


1  Main  menu  11  Selects  gage  calibration  mode,  assumes 

already  logged  on  &  In  main  menu  of 
Process  IDCAL 


2  Cal  range  input  A  A-All  gages 

K-Keyboard  for  individual  gages 
D-Diskfile  for  individual  gages 

3  Cal  range  10  10-Jordan  gages 

3. 25-small  gages 

4  OK  to  dip  Y  Y-Begin  dip  of  rods 

N-Return  to  Step  3 

3  UK  to  cal  0  0-0K  to  wet  gages  for  setting  zero 

1-abort  initial  dip, 

adjust  zero  on  gages  prior  to  activat¬ 
ing  w/rotary  "Balance  R"  on  instrumen¬ 
tation  control  panel,  automatically 
prints  cals  on  LPAO: 

6  Re-do  selected  gages  N  N-No,  usually  OK 

Y-Yes,  if  necessary  to  recal, 

see  Step  6a,  might  wipe  gage  rods  with 

alcohol 

6a  Gage  #'s  to  recal  1-3,  7  Enter  gage  #'s  as  Instructed 

7  Main  menu  15  Exits  process  and  saves  to  file  speci¬ 

fied  as  output 


D4 


Table  D3 

Quadratic  Fit  Calibration  Coefficients  for  Small  Wave  Gages 
Monochromatic  Performance  Teats  in  DSWG  Baaln 


Maximum  Deviation,  10  ft 


Gage  No. 


No. 

Date 

Time 

RD01 

RD02 

RD03 

RD04 

RD05 

RD06 

RD07 

RD08 

RD09 

Min 

Max 

Avg 

1 

18  Feb 

1427 

181 

144 

170 

177 

265 

191 

166 

121 

179 

121 

265 

177 

2 

20  Feb 

1715 

48 

204 

245 

305 

345 

67 

332 

237 

329 

48 

345 

235 

3 

22  Feb 

934 

106 

167 

214 

338 

356 

131 

238 

231 

282 

106 

356 

229 

4 

1700 

68 

84 

143 

295 

254 

102 

203 

130 

206 

68 

295 

165 

5 

23  Feb 

1230 

89 

104 

168 

274 

241 

104 

196 

127 

211 

89 

274 

168 

6 

1600 

81 

122 

169 

265 

244 

144 

217 

147 

228 

81 

265 

180 

7 

25  Feb 

1015 

126 

116 

127 

173 

237 

136 

208 

148 

235 

116 

237 

167 

8 

1900 

99 

79 

134 

185 

234 

115 

195 

146 

255 

79 

255 

160 

9 

27  Feb 

830 

169 

110 

148 

120 

188 

182 

201 

141 

171 

110 

201 

159 

10 

1320 

167 

121 

95 

107 

160 

169 

131 

95 

132 

95 

169 

131 

11 

1900 

78 

87 

142 

150 

159 

111 

175 

127 

184 

78 

184 

135 

12 

4  Mar 

830 

22 

41 

48 

72 

86 

80 

113 

113 

129 

22 

129 

78 

13 

1330 

34 

48 

67 

78 

88 

61 

90 

76 

109 

34 

109 

72 

14 

1930 

37 

51 

56 

72 

65 

51 

85 

81 

114 

37 

114 

68 

15 

6  Mar 

800 

72 

73 

90 

102 

117 

85 

101 

79 

122 

72 

122 

93 

16 

1230 

105 

134 

130 

136 

117 

186 

149 

94 

113 

94 

186 

129 

17 

1920 

20 

35 

50 

55 

56 

45 

54 

62 

63 

20 

62 

49 

18 

8  Mar 

800 

26 

39 

41 

54 

62 

48 

58 

48 

46 

26 

62 

47 

19 

1300 

35 

49 

44 

65 

71 

86 

76 

70 

61 

35 

86 

62 

20 

1600 

87 

97 

81 

49 

74 

51 

59 

40 

58 

40 

97 

66 

21 

9  Mar 

730 

39 

34 

39 

28 

60 

35 

28 

35 

39 

28 

60 

37 

22 

11  Mar 

830 

161 

169 

156 

174 

211 

161 

165 

174 

178 

156 

211 

172 

Averages: 

84 

96 

116 

149 

168 

106 

147 

115 

157 

71 

186 

126 

Note:  DSWG  **  Directional  apectral  wave  generator. 


Table  D4 

Quadratic  Fit  Calibration  Coefficients  for  Jordan  Wave  Gages 
Monochromatic  Performance  Tests  in  DSWG  Basin 
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Note:  DSWG  -  Directional  spectral  wave  generator. 


appendix  E:  WAVE  generation  and  measurement 


Table  El 

Procedure  for  Process  TAPEM2 
Wave  Generation  and  Measurement 


Description 

Input 

Comments 

Log  on 

User  ID: 

Password : 

TAPEM2 

VAX 

Get  current  password  from  Hampton 

Tape  drive 

0 

O-MSAO : ,  1-MSBO: 

Control  signal 
filename 

M17360 

Must  be  6  characters  long  &  exist  on  tape 
Indicated  In  Step  2 

Starting  delay  to 
disk  I/O,  sec 

26 

Time  after  starting  test  before  disk  storage 
of  measured  data  begins.  Includes  10-sec 
delay  time 

Total  test  length, 
sec 

25 

Total  time  DSWG  will  run  and  data  will  be 
stored  on  disk,  make  GT  time  specified  In 
Header  2  of  Process  IDCAL 

Options 

S 

S-Starts  test 

E^Ends  test  if  abort,  enter  positive  real  # 
for  #  of  seconds  to  ramp  down,  normally  ends 
automatically  after  time  specified  in  Step  5 
above 

Other  options  not  used  include: 

D-Change  channels  saved  on  disk  or  add 
decimation 

0«Graph  channels  on  oscilloscope 

P**Print  data  on  Versatec  plotter 

Note:  DSWG  “  Directional  spectral  wave  generator. 
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appendix  F:  wave  analysis 


O  >-  >-  >  >-  N  CM 


Command  File  for  Program  DISRDISK 


*  SET  VERIFY 

*  SDA 

*  R  DISKDISK 
Ml 7360 

1 

*  PU  WVANSDISK.COM 

*  SET  NOVERIFY 


Command  File  for  Program  WVANSDISK1 


*  SET  VERIFY 

*  ON  ERROR  THEN  GOTO  EXIT 

*  R  WVANSDI SKI 
-1 


Ml 7360 ,  D*30,  S*6' 

1.0 


6  @  T=2.00 


Ml 7360 

♦  PRINT  FOR001.DAT 

*  PU  FOR001.DAT 

*  PU  [HAMPTON3DISKDISK.COM 
!*  REN  PARM.PLV  M17360.PPLV 
!*  REN  VECTR1.PLV  M17360.VPLV 

♦  PHASE2 

♦  EXITi 

*  SET  NOVERIFY 


S 


■t>M2ZZZ-<o 


Command  File  for  Program  WAVELS 


C 

*  SET  VERIFY 

*  ON  ERROR  THEN  GOTO  EXIT 

*  R  [TURNER] WAVELS 
-1 


17360, 

.000 


DIR=30,  S=6",  6  <?  T=2 . 00  SEC 


Ml 7360 

*  PRINT  FOR001.DAT 

*  PU  FOR001.DAT 

*  EXIT! 

*  SET  NOVERIFY 


Table  FI 


Procedure  for  Program  DISKDISK 
Stage  1,  Preprocessing  Wave  Analysis 


Step 

Description 

Input 

Comments 

1 

Log  on 

User  ID: 

Password : 

Hampton 

Get  current  password  from  Hampton 

2 

End  of  test 

End 

Sends  message  to  users  informing  testing  Is 
completed 

3 

Preprocessing  edit 

EDD 

Edits  command  file  DISKDISK,  Input  de¬ 
scribed  in  Steps  4  &  5,  close  file  using 
standard  procedure 

4 

Disk  filename 

M17360 

5 

Run  # 

1 

Enter  appropriate  run  # 

6 

Run  DISKDISK 

RDD 

Combines  calibration  coefficients  from 
Process  IDCAL  w/data  from  Process  TAPEM2 
into  dlskfile  specified  in  Step  4.  Now  In 
[DATA  ANAL]  subdirectory 

S£££ 

1 

2 

3 

4 

5 

6 

7 

8 


8a 

8b 

9 

9a 

9b 

10 

11 


Table  F2 

Procedure  for  Program  WVANSDISK1 
Stage  2,  Zero-Crossing  Wave  Analysis 


Description 

Input 

Comments 

Log  on 

User  ID: 

Password : 

Hampton 

Get  current  password  from  Hampton 

Change  directory 

SDA 

Set  default  to  [DATA_ANAL].  If  run  con¬ 
currently  w/Stage  1,  already  in  subdirec¬ 
tory,  omit  step 

Zero-crossing  edit 

EWA 

Edits  command  file  WVANSDISK1.  Input  de¬ 
scribed  in  Steps  4-12 

Analysis  periods 

-1 

-1-All 

Enter  #  periods  .LE.  #  collected  &  speci¬ 
fied  in  Process  IDCAL 

Skip  records 

0 

O-process  all 

Enter  #  records  to  skip  prior  to  processing 

Limited  output 

Y 

Y-Std  output,  N-Detailed  output 

Summary  output 

Y 

Y-Std  output,  N-Detailed  output 

Process  channels 

Y 

Y-All,  skip  Steps  8a-8b  below 

N-Process  selected  channels  only,  see 

Steps  8a-8b 

How  many 

3 

Enter  #  channels  to  process,  .LE.  total 
available 

Channel  I's 

O  00  M 

Enter  channel  #'s  to  process,  one  per  line 

Plot  data 

Y 

Y-Yes,  N-No  plots  desired,  skip 

Steps  9a-9b 

x-axls  length 

7. 

Enter  length  of  x-axls  in  inches 

x-axls  increment 

2. 

Enter  x-axis  increment  in  seconds 

Descriptive  title 

Enter  .LE.  40  characters 

Depth 

1.0 

Feet,  does  not  seem  to  make  any  difference 

(Continued) 
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Table  F2  (Concluded) 


Step 

Description 

Input 

Comments 

12 

Filename 

Ml  7360 

Enter  dlskfile  containing  data,  same  as 
created  In  Stage  1 

13 

Run  WVANSDISK1 

RWA 

Runs  WVANSDISK1,  zero-crossing  average  & 
significant  wave  periods  &  heights  for 
model  &  prototype.  Prints  tabular  lists 
plots  (if  requested) 

Table  F3 

Procedure  for  Program  WAVELS 
Stage  3,  Harmonic  Wave  Analysis 


Step 

Description 

Input 

Comments 

1 

Log  on 

User  ID: 

Password : 

Hampton 

Get  current  password  from  Hampton 

2 

Change  directory 

SDA 

Set  default  to  [DATA  ANAL] 

3 

Harmonic  analysis 

EHA 

Edits  command  file  WAVELS. 

Input  described  in  Steps  4-13 

4 

Analysis  periods 

-1 

-  1-All 

Enter  #  periods  .LE.  if  collected  &  speci¬ 
fied  in  Process  IDCAL 

5 

Skip  records 

0 

0-process  all 

Enter  #  records  to  skip  prior  to 
processing 

6 

Process  channels 

Y 

Y-All,  skip  Steps  6a-6b  below 

N=Process  selected  channels  only,  see 

Steps  6a-6b 

6a 

How  many 

3 

Enter  #  channels  to  process,  .LE.  total 
available 

6b 

Channel  #'s 

7 

8 
q 

Enter  channel  if' s  to  process,  one  per  line 

7 

Plot  raw  data 

7 

N 

Y-Yes,  N-No  plots  desired 

8 

Plot  raw  data  &  fit 

N 

Y-Plots  of  raw  &  least-squares  fit  data 
desired 

N-No  plots  desired 

9 

Plot  residuals 

N 

Y-Plots  of  residuals  desired 

N-No  plots  desired 

10 

Descriptive  title 

Enter  .LE.  40  characters 

11 

Wave  period 

2.0 

Fundamental  period  in  seconds 

12 

#  Components 

4 

Enter  #  of  harmonics  desired 

(Continued) 


Table  F3  (Concluded) 


Step 

Description 

Input 

Comments 

13 

Disk  filename 

Ml  7360 

Enter  diskfile  containing  data,  same  as 
created  in  Stage  1 

14 

Run  WAVELS 

RHA 

Runs  [TURNER] WAVELS  harmonic  analysis  for 
percent  variance  in  each  harmonic. 

Prints  tabular  lists  &  plots  (if 
requested) 

15 

Logoff 

LO 

F9 


o 

B 

*3 

c 


g 

E 

f 

f  (w) 

Vf) 

F3(f,e) 

h 

H 

H 

«b 

H 

n 

H 

P 

H/L 

H/S 

H(0) 


J 

J 

k 

L 


Real  Fourier  coefficients 
Amplitude  of  component 
Fourier  coefficient 
Mean  water  surface  elevation 
Width  of  paddles  of  DSWG 
Imaginary  Fourier  coefficients 
Wave  celerity  or  phase  speed 
Group  velocity  of  wave 
Variance  of  wave  signal 
Frequency,  Hz 

Function  of  0  with  w  as  a  dummy  variable 
Two-dimensional  wave  height  transfer  function 
Three-dimensional  wave  height  transfer  function 
Water  depth 
Wave  height 

Average  wave  height  in  measurement  area 

Maximum  prebreaking  wave  height 

Wave  height  measured  at  gage  n 

Wave  height  of  the  pC^  spurious  wave  component 

Wave  steepness 

Helght-to-stroke  ratio 

Wave  height  of  main  component  wave;  wave  height  corrected  for 
directional  effects 

.th 

j  component 

Total  number  of  components 
Wave  number 

Linear,  shallow-water  wavelength 


G3 


N 

N! 

P 

S 

sc(y.O 

t 

T 

T 

min 

T1 

T2 

v 

X 

Y 

ah 

AL 

At 

Ay 

£<t) 

e(t) 

n(t) 

n  (t) 
in 

njt) 

0 

0 

m 

0 

P 

n 


Total  number  of  data  samples 

Number  of  DSWG  paddles  required  to  make  one  cycle 
Index  for  spurious  wave  components 
Wavemaker  stroke 

Control  signal  for  linear  sinusoidal  waveforms 
Time 

Time  for  wave  to  travel  between  two  points 

Minimum  period  below  which  spurious  waves  will  be  generated 

Travel  time  to  back  gage  row  located  30  ft  from  DSWG 

Travel  time  to  reflect  off  back  wall  and  return  to  back  gage  row 

Dummy  variable  of  Integration 

Coordinate  axis 

Distance  along  DSWG  corresponding  to  one  cycle;  coordinate  axis 
Mean  wave  height  variation 
Incremental  wavelength 
Time  interval 

Distance  between  wave  gages  parallel  to  DSWG 
Paddle  displacement  time  series 

Noise  component  of  water  surface  elevation  time  series 

Water  surface  elevation  time  series 

Measured  water  surface  elevation  time  series 

True  component  water  surface  elevation  time  series 

Wave  direction,  angle  of  wave  propagation;  angle  between  wave  crest 
and  DSWG 

Measured  wave  direction 

Direction  of  pC^  spurious  wave  component 

Goda’s  nonlinear  parameter 


G4 


3.14159 


Phase  lag  between  paddles  of  DSVG 
Phase  of  j*'*1  component 

Total  phase  lag  tor  wave  generation,  consists  of  frequency  increment, 
constant,  and  offset  portions 

Constant  plus  offset  phase  lag  for  wave  generation 
Angular  frequency  of  J**1  wave  component 


